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ABSTRACT 


A  mathematical  model  is  developed  to  analyze  the  seismic 
response  of  a  drydocked  vessel  in  three  degrees  of  freedom; 
vessel  rotation  about  its  keel  and  vessel  horizontal  and 
vertical  translations  relative  to  the  drydock  cradle.  Data 
from  eleven  actual  vessel -drydock  systems  and  the  time 
acceleration  history  of  an  earthquake  are  implemented  to 
predict  vessel  three  degree  of  freedom  response  during  an 
earthquake.  Vessel  seismic  response  and  the  resultant  drydock 
forces  are  compared  to  vessel-drydock  system  failure  criteria 
to  determine  stability  of  the  system  during  earthquakes.  The 
three  degree  of  freedom  vessel  response  model  is  compared  to  a 
one  and  two  degree  of  freedom  vessel  response  models  and  a 
model  in  which  seismic  loading  is  simulated  by  a  single  static 
force.  The  three  degree  of  freedom  vessel  motion  is  shown  to 
be  the  most  accurate  method  for  analyzing  vessel-drydock 
system  failure  criteria., 
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1 . 0  INTRODUCTION 


1.1.  Background 

In  the  design  of  a  drydock  for  military  vessels,  an 
important  factor  is  the  environmental-generated  forces  that 
are  found  in  the  drydock  cradle  which  supports  the  ship. 
Earthquakes  and  the  wind  are  considered  the  two  most 
significant  environmental  loading  factors.  This  thesis 
explores  the  seismic  loading  aspect  of  the  problem. 

Regardless  of  its  location,  any  U.S.  Navy  ship  that  is 
nuclear  powered  and  in  drydock  is  considered  to  be  in  a  high 
seismic  risk  area.  Currently  the  Navy  defines  a  high  seismic 
risk  area  as  one  in  which  earthquake  forces  be  approximated  by 
a  steady  horizontal  force  of  0.2  g  times  the  vessel  mass 
acting  at  the  center  of  gravity.  The  drydock  can  be  analyzed 
to  see  if  it  can  withstand  this  quasi-static  loading. 

1.2  Previous  Work 

In  1981  B.V.  Viscomi  studied  the  seismic  response  of  a 
drydocked  submarine  using  a  single  degree  of  freedom  model 
[14].  This  analysis  assumed  that  the  vessel  was  allowed  to 
rotate  about  the  keel.  Viscomi  used  the  time  history 
acceleration  record  of  the  North-South  component  of  the  1946 
El  Centro,  CA  Earthquake,  California  Institute  of  Technology 
processing  scheme.  The  vessel  was  analyzed  for  system  failure 
(i.e.  the  vessel  lifting  off  of  a  row  of  side  blocks). 

A  thesis  by  C.F.  Barker  in  1985  [1]  used  a  two  degree  of 
freedom  model  to  study  seismic  response.  Employing  the  1946 
El  Centro  Earthquake,  Massachusetts  Institute  of  Technology 


standard  processing  scheme,  Barker  analyzed  rotation  about  the 
keel  plus  horizontal  relative  translation  between  the  vessel 
and  the  drydock  cradle  perpendicular  to  the  vessel  longitudal 
axis.  The  two  degree  of  freedom  model  permitted  the  analysis 
of  drydock  cradle  failure:  block  sliding  and  tipping,  which 
were  not  possible  to  analyze  in  Viscomi's  one  degree  of 
freedom  model,  along  with  block  liftoff. 

1.3  Contributions  of  This  Thesis 

This  thesis  develops  a  three  degree  of  freedom  model  of  a 

vessel  in  drydock.  Motion  will  be  permitted  in  three 
directions:  rotation  about  the  keel  and  relative 
translations,  horizontal  and  vertical,  between  the  vessel  and 
the  drydock  cradle.  The  three  degree  of  freedom  model  will 
make  possible  the  analysis  of  drydock  cradle  failure: 
additional  block  crushing  forces  due  to  vertical  motion,  which 
was  not  conceivable  in  the  previous  one  and  two  degree  of 
freedom  models  developed  by  Viscomi  and  Barker  respectively, 
along  with  block  sliding,  tipping  and  liftoff. 

The  seismic  input  for  the  three  degree  of  freedom  model 
is  the  time  acceleration  history  of  the  North-South  component 
of  the  1946  El  Centro,  CA  Earthquake,  Massachusetts  Institute 
of  Technology  standard  processing  scheme.  The  MIT  version  of 
the  El  Centro  Earthquake  acceleration  record  is  described  in 
reference  [9]  ,  and  is  displayed  graphically  in  Figure  1.1. 

1.4  Outline  of  This  Thesis 

Section  2  of  this  thesis  describes  the  vessel-drydock 


system  for  the  reader  and  introduces  terms  which  are  used 
throughout  the  thesis.  Section  3  describes  the  various 


failure  mechanisms  found  in  the  vessel-drydock  system  which 
can  occur. 

Section  4  examines  the  modelling  of  the  vessel-drydock 
system  under  seismic  loading.  Section  4.1  describes  the 
quasi-static  loading  method,  and  Section  4.2  develops  the  one 
degree  of  freedom  equation  of  vessel  motion  model.  Section 
4 . 3  develops  the  two  degree  of  freedom  equations  of  vessel 
motion  model.  Finally,  Section  4.4  develops  the  three  degree 
of  freedom  equations  of  motion  model. 

Section  5  discusses  the  system  parameters  that  are 
required  for  the  mathematical  model  found  in  Section  4 . 
Included  in  this  Section  is  the  modelling  of  block  stiffness. 
Eleven  typical  vessel-drydock  system  configurations  are 
studied  to  be  implemented  into  Section  6. 

Section  6  examines  methods  of  evaluating  the  one,  two  and 
three  degree  of  freedom  models  found  in  Section  4.  Section 
6.1  evaluates  the  linear  equations  of  motion  models  using 
modal  analysis  method.  Section  6.2  evaluates  both  linear  and 
non-linear  equations  of  motion  models  using  a  fourth-order 
Runge-Kutta  numerical  analysis  scheme.  Section  6.3  explains 
the  response  spectrum  method  of  determining  the  maximum  value 
of  vessel  seismic  response.  Section  6.4  determines  system 
response  using  quasi-static  force  method. 

Section  7  describes  the  development  and  testing  of  the 
one,  two  and  three  degree  of  freedom  vessel  motion  computer 
programs,  and  contains  the  results  for  a  seismic  analysis  of 
several  actual  vessel-drydock  systems. 

Section  8  summarizes  the  response  predictions  for  several 
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actual  vessel-drydock  systems  obtained  with  quasi-static,  one 
degree  of  freedom,  two  degree  of  freedom  and  three  degree  of 
freedom  analysis  methods,  and  conclusions  are  drawn. 

Appendix  1  is  an  example  listing  of  the  fourth-order 
Runge-Kutta  computer  program.  Appendix  2  is  an  example 
listing  of  the  modal  analysis  computer  program.  Finally, 
Appendix  3  describes  the  modal  analysis  of  the  two  and  three 
degree  of  freedom  vessel-drydock  systems  and  predicts  maximum 
system  response  using  the  response  spectrum  method  of  analysis 
and  modal  participation  factors. 


ACCELERATION  IN  CM/SEC 


In  order  to  perform  inspection,  maintentance  and  repair 
on  the  outer  hull  of  a  ship,  it  is  necessary  to  place  the 
vessel  in  a  drydock.  A  drydock  is  a  concrete  encasement  built 
into  the  earth.  At  one  end  of  the  drydock  is  a  flood  gate  to 
allow  the  entry  and  departure  of  the  vessel.  A  cradle,  formed 
by  drydock  blocks,  is  built  upon  the  drydock  floor  in  order  to 
support  the  ship  once  the  water  is  pumped  out  of  the 
encasement.  The  drydock  must  be  able  to  hold  the  weight  of 
the  vessel  and  cradle. 

Due  to  the  various  hull  configurations  associated  with 
the  different  classes  of  ships,  the  cradle  component  of  the 
drydock-cradle-vessel  system  is  varied  to  suit  a  particular 
docking  situation.  The  cradle  accomplishes  the  transfer  of 
vessel  weight  to  the  drydock  floor,  provides  stable  support  to 
the  vessel  and  allows  access  to  the  vessel  hull.  This  cradle 


is  usually  constructed  of  timber,  concrete,  or  a  timber- 
concrete  composite. 


System  Component 


The  primary  components  of  the  drydock-cradle-vessel 


system  are  shown  in  Figure  2.1.  This  section  defines  terms, 
related  to  docking,  used  in  this  thesis  to  describe  the  system 
components . 


Blocks  -  The  units,  consisting  of  timber,  concrete,  steel 
and  other  materials,  which  together  make  up  the  cradle. 

Cradle  -  A  framework  of  blocks  which  supports  the  vessel 
when  the  drydock  is  dewatered. 

KG  -  The  distance  between  the  vessel  baseline  (ship’s 
keel)  to  its  vertical  center  of  gravity. 

Keel  Blocks  -  The  center  blocks,  directly  beneath  the 
vessel's  keel. 

Pier  -  A  column  built  of  blocks  that  extends  from  the 
ship  hull  to  dock  floor. 

Side  Blocks  -  The  blocks  located  to  the  right  and  left  of 
the  keel  blocks. 

Ton  -  A  long  ton,  2240  lbs. 

g-  Acceleration  of  gravity,  32.2  ft/sec2. 


3.0  VESSEL-DRYDOCK  SYSTEM  FAILirRF  HECHANISMS 


There  are  four  failure  modes  that  the  vessel-drydock 
system  can  exhibit:  block  crushing,  block  sliding,  block 
overturning  and  block  liftoff  (i.e.  the  liftoff  of  vessel  from 
a  row  of  blocks) . 

3.1  Block  Crushing6 

A  block  can  support  compressive  stresses  linearly  until 
the  blocks'  proportional  limit  is  exceeded.  The  proportional 
limit  is  a  material  property  of  the  block  and  is  the  maximum 
compressive  stress  at  which  stress  is  still  linearly 
proportional  to  strain.  If  the  compressive  proportional  limit 
is  exceeded,  the  block  is  considered  to  have  failed. 

3.2  Block  Sliding6 

Due  to  friction  at  the  interfaces  between  the  blocks  and 
the  vessel  and  between  blocks  with  other  blocks,  drydock 
blocks  will  have  the  tendency  to  resist  sliding  when  subjected 
to  vertical  and  horizontal  loads.  This  is  true  when  in  the 
absence  of  mechanical  fasteners. 

First  consider  the  blocks  that  form  the  keel  pier. 

Figure  3.1  shows  that  the  resistance  of  the  keel  pier  to 
sliding  is 

H  =  UjV  (block-block  interface) 

and  H  =  P2V  (vessel-block  interface) 


where 


H  =  horizontal  resisting  force 
y -j_  =  coefficient  of  friction,  block-block  interface 
y2  =  coefficient  of  friction,  block-vessel  interface 
V  =  vertical  load  on  the  blocks. 
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Failure  of  the  keel  pier  due  to  sliding  will  occur  when 
the  ratio  of  horizontal  force  to  the  vertical  force  is  greater 
than  the  lowest  coefficient  of  friction  found  in  the  block.  In 
other  words,  sliding  occurs  when 
|  >(Vi  or  y2)  • 

As  with  the  keel  piers,  the  resistance  of  the  side  piers 
to  sliding  will  depend  on  the  interface  which  offers  the  least 
resistance,  the  block-vessel  interface  or  the  block-block 
interface.  However,  this  case  is  more  complicated  due  to  the 
complex  geometry  of  the  side  pier  as  shown  in  Figure  3.2. 

The  block-vessei  interface  of  the  side  pier  is  examined 
and  shown  in  Figure  3.3.  An  arbitrary  force  F  at  some  angle 
from  vertical  is  applied  to  the  face  for  the  side  block  cap. 

At  this  interface,  the  horizontal  and  vertical  block  force 
reactions  are: 

Normal  Force  =  H  sin®  +  V  cos <J> 

Tangential  Force  =  H  cos®  -  V  sin® 
where  H  =  Horizontal  Reaction 

V  =  Vertical  Reaction 
®  =  Block  Inclination  Angle. 

Using  the  formula  for  horizontal  sliding  resistance,  H  =v 
V,  previously  introduced  in  the  keel  pier  case  and  the 
above  equations  yields 

H  cos®'-  V  sin®  =  y2  (H  +  v  cos®)  . 

Rearranging  the  above  equation  gives 

H  _  vi 5  cos®  +  sin® 


where 


y2  cos4>  +  sin4> 

^2  ~  cos$  -  y2sin^ 

For  the  block-block  interface  of  the  side  pier,  sliding 
resistance  is  found  by  taking  the  vertical  and  horizontal 
components  of  the  applied  force,  F,  and  using  them  in  the 
formula  H  =  y  x  V.  Like  in  the  keel  pier,  failure  of  the  side 
pier  due  to  sliding  will  occur  when  the  following  is  satisfied 
~  >  (Pi  or  P2  ' )  • 

There  is  a  critical  angle  a  where  no  slippage  occurs. 
Referring  again  to  Figure  3.3,  this  angle  can  be  calculated  by 
comparing  the  maximum  tangential  force  F^  max  the  interface 
can  have  without  slipping  and  the  corresponding  normal  force 
Fn.  The  following  relationship  holds 
Ft  max  =  ^Fn 


Finally, 

a  =  arc  tan  y. 

Since  the  applied  force  F  acts  in  a  straight  line  through  all 
blocks  to  the  ground,  the  angle  of  for  each  interface  must  be 
calculated.  If  force  F  is  applied  outside  of  these  angles, 
slippage  will  occur. 

3.3  Block  Overturning11 

The  third  failure  mechanism  is  the  overturning  of  a  block 
due  to  an  applied  force.  The  line  of  action  of  this  force 
must  fall  within  the  middle  one-third  of  the  base  of  the 
block,  as  shown  in  Figure  3.4,  or  the  block  will  tip  over. 

In  order  to  show  the  region  that  the  block  will  remain 
upright,  the  inclined  force  F  must  be  broken  into  a  transverse 


component  and  a  vertical  component  Fv.  Superposing  the 
bending  and  axial  stresses  caused  by  F^  and  Fv  gives  the 
resultant  stress  at  any  point  in  the  block  as 


where  L  =  block  height 

A  =  base  area,  here  bh 
Iz  =  base  moment  of  intertia,  here  bh3/l2. 

The  minimum  stress  point  on  the  bottom  of  the  block 
occurs  at  point  M  (Figure  3.4),  a  distance  -h/2  from  the  z 
axis.  Hence, 

F 

In  \  =  Ik.  P  _ 

°x  m  bh  h  bh  * 

When  °xitl  at  M  is  negative,  compressive  stress  is  present  and 
the  block  is  stable.  When  at  M  is  positive,  tensile 
stress  is  present  and  the  block  will  overturn  since  no 
fasteners  hold  block  to  the  ground  to  develop  tensile  stress. 
The  limit  of  block  stability  is  when 

(°x)m  =  0> 

Thus,  the  following  equation  holds 

F  =  —  F 
h  6L  v* 

This  condition  exists  when  force  F  is  applied  at  an  angle  such 
that  the  line  of  force  lies  within  the  one-third  of  the  base 
of  the  block.  If  the  line  of  force  is  outside  this  region, 
the  block  will  tip. 

3.4  Vessel  Liftoff14 

The  fourth  and  final  failure  mode  occurs  when  the  vessel 
breaks  contact  with  either  the  side  or  keel  piers.  This 
failure  occurs  when  the  dynamic  deflection  of  a  row  of  blocks 
is  equal  to,  or  exceeds  the  average  static  deflection  of  the 


blocks.  The  static  deflection  of  blocks  is  given  by 


where 


6s  = 


w 


2K 


sv 


+  K 


kv 


W  =  submarine  weight 
Ksv  =  side  pier  vertical  stiffness 
Kkv  =  keel  pier  vertical  stiffness, 


The  dynamic  deflection  can  be  due  to  vertical  displacement  y, 
rotational  displacement  6,  or  a  combination  of  both. 


vs 
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4.0  MATHEMATICAL  MODELS  OF  VESSEL-DRYDOCK  SYSTEM 

This  portion  of  the  thesis  discusses  the  modelling  of  the 
drydocked  submarine  subjected  to  a  seismic  load.  Models  will 
include  a  quasi-static  method,  one,  two  and  three  degree  of 
freedom  mathematical  models.  The  models  developed  in  this 
section  will  be  analyzed  in  Section  6  and  7  to  predict  the 
vessel  seismic  response. 

4.1  The  Quasi-static  Method  of  Modelling  Seismic  Response 

4.1.1  Approach5 > 6 

Current  U.S.  Navy  design  method  for  submarine's  seismic 
response  will  be  examined  in  this  section.  The  submarine  is 
described  as  a  rigid  cylindrical  body  with  its  weight  evenly 
dispersed  longitudinally. 

The  quasi-static  force  method  replaces  the  earthquake 
motions  by  a  force  corresponding  to  the  vessel  mass  time 
0.2  g.  This  force  is  horizontally  applied  to  the  submarine's 
center  of  gravity  in  the  tranverse  direction.  The  drydock 
blocking  system  is  determined  by  the  evaluation  of  the  loads 
generated  by  the  quasi-static  force. 

4.1.2  Force  Equations 

The  application  of  the  static  force  to  the  center  of 
gravity  of  the  submarine  in  drydock  is  represented  in  Figure 
4.1.  The  seismic  overturning  moment,  Ms,  is  defined  as 

Ms  =  (A/g) (a) (KG) (2240)  ft-lb 

where 


A  =  vessel  displacement  in  long  tons 

g  =  acceleration  of  gravity 

a  =  vessel's  center  of  gravity  acceleration 


Current  U.S.  Navy  design  practice  states  that  the  acceleration 
of  the  submarine's  center  of  gravity,  a,  should  be  set  equal 
to  0.2  g.  Hence 

Ms  =  448(A) (KG)  ft-lb 

After  determining  the  seismic  overturning  moment,  the 

number  of  sideblocks  required  on  one  side  of  the  vessel  is 

M 

s 


where  N  =  #  of  side  blocks  required  in  the  row 

A  =  contact  area  of  a  block  in  square  inches 

Sp  =  proportional  limit  of  the  block  cap  in  lb/in2 

L  =  distance  between  centerline  of  keel  and  side 
blocks,  in  feet,  as  shown  in  Figure  4.1. 

4.1.3  Vessel  Response  in  the  Quasi-static  Load 
Method 

The  application  of  a  force  applied  at  a  vessel's  center 
of  gravity  will  generate  reactionary  side  block  forces  in  the 
vertical  direction  to  oppose  the  seismic  moment,  as  shown  in 
Figure  4.2.  Summation  of  moments  about  the  keel  yields  the 
equation 

IM  =  Ms  -  L(Fsr  +  Fmr)  +  L(Fsl  -  Fml)  =0  (4.1.1) 


where 

Fs  =  static  pier  forces  (left  and  right) 

Fm  =  pier  forces  due  to  the  applied  moment  (left  and  right), 
Due  to  symmetry  in  the  vessel-drydock  system,  the 


applied  moment  resistive  forces  in  the  side  blocks  are  equal 


in  magnitude  but  opposite  in  direction.  Thus,  Equation  4.1.1 
simplifies  to  the  form 


Ms  =  2L(Fm)  (4.1.2) 

where  Fm  =  Fml  =  Fmr. 

The  force  Fm  is  equal  to  the  side  pier  compression 
displacement,  <5  m,  times  its  vertical  stiffness,  Kvs.  This 


compression  displacement  can  be  expressed  as 
=  L  sin (6)  . 

where  0  =  angle  or  vessel  rotation  about  its  keel. 


(4.1.3) 


Assuming  the  rotational  angle  is  small,  Equations  4.2.1  and 
4.1.3  are  combined  to  yield 


2 

2L  K 


^  *  4 


(4.1.4) 


4.2  One  Degree  of  Freedom;  Rotational  Response 
4.2.1  Approach 

A  one  degree  of  freedom  mathematical  model  of  a  submarine 
in  drydock  is  developed  in  this  section.  The  submarine  is 
identified  as  a  rigid  cylindrical  body  with  an  even 
longitudinal  weight  distribution.  The  system  damping  and 
stiffness  are  provided  by  the  drydock  blocking  arrangement. 
System  excitation  is  due  to  seismic  ground  accelerations. 

Other  assumptions  are 


-» 


-  The  keel  and  side  piers  remain  vertical  during  the 
ground  motions, 

-  no  slippage  between  the  block  cradle  and  drydock,  and 

-  the  keel  and  side  piers  bases  accelerate  at  the  same 
rate  as  the  drydock  (ground) . 

The  idealized  model  of  the  vessel-drydock  system  is  shown 
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in  Figure  4.3.  In  this  system,  vessel  motion  is  restricted  to 
rotation  about  the  keel.  The  side  piers  are  modelled  as 
vertical  due  to  typically  small  angles  or  inclination  of  the 
tops  of  the  side  piers.  Only  vertical  displacements  in  the 
side  piers  are  considered. 

The  one  degree  of  freedom  mathematical  model  predicts  and 
analyzes  vessel  rotational  response  about  its  keel  caused  by 
seismic  ground  accelerations  normal  to  the  submarine's 
longitudinal  axis.  The  resultant  angle  or  rotation,  expressed 
in  radians,  and  side  pier  forces,  generated  by  the  rotation, 
are  compared  to  applicable  pier  failure  criteria  found  in 
Section  3,  to  see  if  the  blocking  arrangment  will  remain 
intacted. 


4.2.2  One  Degree  of  Freedom  Equation  of  Motion  Model 
For  the  one  degree  of  freedom  case,  the  keel  is  the 
origin  with  clockwise  rotation,  0,  positive.  The  simplified 
vessel-block  system  showing  coordinate  0  at  rest  and  when 
excited  are  shown  in  Figures  4.3  and  4.4  respectively.  The 
moment  equilibrium  equation  about  the  keel  (origin)  is: 


^Mk  =  Ik0  +  MKGXg  -  MKGyg  sin© 


=  B/2 (Flsv  -  Fldv)  -  B/2(F2sv  +  F2dv)  +  WKG  sin e  (4.2.1) 


where  Ik  =  Mass  moment  of  inertia  of  vessel  about  the  keel 
M  =  Vessel  Mass 

Flsv'F2sv  =  static  side  block  forces  in  the  vertical  direction 
(Figure  4.4) 

Fldv'F2dv  =  dynamic  side  block  forces  in  the  vertical 


direction  (Figure  4.4) 

W  =  vessel  weight  (Mg) 

and  Xg,yg  =  seismic  ground  accelerations  in  their  respective 
directions 

This  equation  can  be  simplified  by  the  following: 

1. )  The  left  and  right  side  block  rows  have  the  same 

vertical  spring  constant,  Ksv,  have  the  same  static 
deflection  and  equal  but  opposite  dynamic  deflection  due 
to  symmetry ,  i . e . 

^Isv  =  ^2sv  ~  ^sv  (4.2.2) 

and  I  F idv I  =  I  B2dvl  =  ^dv  (4.2.3) 

2. )  The  dynamic  side  block  force,  Fdv,  as  seen  from  Figure 

4.5  equals  to  a  force  due  to  the  modelled  spring 
displacement  plus  a  dissipative  force  due  to  system 
damping.  In  other  words, 

*dv(sPring)  =  (B/2)KSv  sin  6 
and  Fdv (dissipative)  =  C(B/2)  6 
or  when  combined, 

Fdv  =  B/2  (Ksv  sin?  +  C 6)  (4.2.4) 

where  Ksv  =  vertical  spring  constant  of  the  side  blocks 

C  =*  vertical  damping  coefficient  of  the  side  blocks 
Substituting  Equations  4.2.2,  4.2.3,  and  4.2.4  into 
Equation  4.2.1  and  rearranging  components  yields 
Ike  +  (B2/2)C6  +  [(B2/2)Ksv  -  W  KG]  sine 

=  -MKG  xg  +  MKG  yg  sin  6  (4.2.5) 

An  additional  simplification  is  that  since  0  is  very  small, 
sin  0  =  6. 

The  final  simplification  is  that  damping  is  expressed  as  a 


fraction  of  critical  damping.  To  do  this,  allow 
(B2/2)C  =  C  Ccr  =  2K  1]^ 
where  £  =  viscous  damping  factor(B2C/2Ccr) 

Ccr  =  critical  damping  coefficient  =  21^'^ 
wn  =  system  undamped  natural  frequency. 

Now,  Equation  4.2.5  simplified  to 

Iks  +  2  ?  IKwne  +  [  (B2/2 )  Ksv  -  WKG]6 

=  -MKG  xg  +  MKG  yg  0  (4.2.6) 

The  non-linear  term,  MKG  yg e,  found  in  the  equation  of 
motion  will  be  removed  at  this  time  to  linearize  Equation 
4.2.6  since  6  is  assumed  to  be  small.  The  effects  of  the 
non-linear  term  will  be  evaluated  in  Section  6  and  7. 

The  final  form  of  the  linearized  one  degree  of  freedom 

equation  of  motion  is 

•  •  • 

mi6  +  Cj6  +  k 2®  =  _ni2  *g  (4.2.7) 

where  =  Ik 

C1  “  2^  ik^n 

kx  =  (B2/2)Ksv  -  W  KG 

m2  =  M  KG. 

In  matrix  form. 


To  simplify  the  explanation,  the  matrices  will  be  redefined  as 
[A] (y )  +  [B] (y }  =  (E(t) }  (4.2.8) 


where 


This  section  develops  a  two  degree  of  freedom 


mathematical  model  of  a  drydocked  submarine.  As  in  the  case 
of  the  one  degree  of  freedom  (Section  4.2.1),  the  model 
assumptions  still  hold,  i.e.  a  rigid  cylindrical  body,  keel 
and  side  piers  remain  vertical,  no  slippage  between  the  block 
cradle  and  drydock,  and  the  bases  of  the  keel  and  side  piers 
accelerate  at  the  same  rate  as  the  drydock  (ground) . 

The  idealized  model  of  the  vessel -drydock  system  is  shown 
in  Figure  4.6.  In  this  system,  vessel  motion  is  restricted 
to  rotation  about  the  keel,  and  in  horizontal  transverse 
translation  relative  to  the  keel  and  side  pier  supports.  The 
approach  for  determining  vessel-drydock  failue  is  identical  to 
the  one  degree  of  freedom  case  except  for  the  addition  of 
sliding  and  tipping  failure  modes. 

The  two  degree  of  freedom  mathematical  model  predicts  and 
analyzes  submarine  rotational  and  transitional  response  caused 
by  seismic  ground  accelerations  normal  to  the  vessel's 
longitudinal  axis.  The  rotational  response  is  in  radians,  the 
translational  response  in  inches,  and  both  responses  are 
expressed  as  a  function  of  time.  The  keel  and  side  pier 
forces  are  compared  to  applicable  pier  failure  criteria  found 


in  Section  3. 


4.3.2  Two  Degree  of  Freedom  Equation  of  Motion  Model 

As  in  the  one  degree  of  freedom  case,  the  keel  is  the 
origin  for  the  two  degree  of  freedom  system  with  clockwise 
rotation,  fc<  ,  is  positive.  The  linear  relative  translation 
coordinate,  x,  is  defined  with  respect  to  vessel  and  drydock 
translation  relative  to  the  ground.  Defining  u  as  the 
position  of  the  submarine  keel  relative  to  the  ground,  then 
the  following  relationships  hold: 

X  =  U  -  Xg 
U  =  X  +  Xg 

and  u=x+Xg.  (4.3.1) 

The  simplified  vessel-block  system  showing  coordinates  (6  and 
x)  at  rest  and  when  excited  are  shown  in  Figures  4.6  and  4.7 
respectively. 

A  balance  of  forces  in  the  x  direction,  as  shown  in 
Figure  4.7,  yields  the  first  equation  of  motion 

IFX  =  Mu  +  MKG  6  =  (Flsh  -  Fldh)  -  (F2sh  +  F2dh) 

~  (F3sh  +  F3dh>  (4.3.2) 

where  F^g^,  F2sfc,  F3sh  =  static  block  forces  in  the  horizontal 

direction  at  their  respective 
position 

Fldh»  F2dh'  F3dh  =  dynamic  block  forces  in  the 

horizontal  direction  at  their 
respective  position. 

The  above  equation  can  be  simplified  by  the  following: 

1)  The  left  and  right  side  block  rows  have  the  same 
horizontal  spring  constant,  Kg^,  have  equal  but  opposite 
static  horizontal  deflection  and  the  same  dynamic  horizontal 


deflection  due  to  symmetry,  i.e., 

Flsh  +  F2sh  =  0 

Fldh  =  F2dh  •  (4.3.3) 

Also,  the  keel  block  rows  do  not  experience  any  net  static 
horizontal  deflection  since  the  inclination  angle  of  the  block 
cap  is  zero.  Hence 

F3sh  =  0-  (4.3.4) 

2)  The  dynamic  block  forces,  F^  as  seen  from  Figure  4.7 
equals  to  a  force  due  to  the  modelled  springs  plus  a 
dissipative  force  due  to  system  damping.  In  other  words 

Fdh  (spring)  =  Fldh  (spring)  +  F2dh (spring)  +  F3dh  (spring) 
=  Kshx  +  Ksh  x  +  Kkhx 
Fdh  (spring)  =  (2Ksh  +  Kkh)x 
Likewise 

Fdh  (dissipative)  =  cxx 
so 

Fdh  =  (2Kj^s  +  Kj^j^)  X  +  CjjX  (4.3.5) 

where  KSh  =  side  pier  stiffness 
Kkh  =  ^eel  pier  stiffness 

Cx  =  system  horizontal  damping  coefficient. 
Substituting  Equations  4.3.1,  4.3.3,  4.3.4,  and  4.3.5 
into  Equation  4.3.2  and  rearranging  components  yields 

mu  x  +  m^2  6  +  ci  x  +  kn  x  =  -mu  Xg  (4.3.6) 

where  m^i  =  M 

®12  =  ^  KG 

Ci  =  Cx  =  system  horizontal  damping  coefficient 
kn  =  (2Ksft  +  Kj^)  =  system  horizontal  stiffness 
Summing  the  moments  about  the  origin  in  a  similar  fashion 


as  the  previous  degree  of  freedom  case,  as  shown  in  Figure 
4.7,  yields  the  second  equation  of  motion, 


IMk  =  Ik6  +  MKG  u  -  MKG  yg  sine 

=  B/2 (Flsv  -  Fldv)  -  B/2(F2sv  +  F2dv)  +  WKG  sin6  (4.3.7) 

The  second  equation  can  be  reduced  as  it  was  in  the  one  degree 
of  freedom  case  to 

Ike  +  MKG  X  +  Ce6  +  (B2/2  Ksv  -  WKG) 6 

=  MKG  Xg  +  MKG  yg6.  (4.3.8) 

Once  again,  the  non-linear  term,  MKG  yg  6,  found  in 
Equation  4.3.8  will  be  removed  leaving  a  linearized  equation 
since  P  is  assumed  to  be  small.  The  effects  of  the  nonlinear 
term  will  be  evaluated  in  Sections  6  and  7. 


Now  the  system  of  equations  for  the  two  degree  of  freedom 
vessel -drydock  system  is  as  follows, 

mi!  x  +  m12t  +  clx  x  +  )cn  x  =  -mlx  xg  (4.3.9a) 

and 

m22e  +  m21  x  +  c12*  +  k22e  =  -m21  xg  (4,3.9b) 

where  =  M 

m2i  —  m12  =  MKG 
®22  =  xk 

c^  =  system  horizontal  damping  coefficient 
c2  =  system  rotational  damping  coefficient 
kn  =  2Ksj1  +  K]^ 
k22  =  B2/2  Ksv  -  WKG. 

These  two  equations  of  motion  are  coupled  in  the  mass  times 
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acceleration  term,  which  is  known  as  inertial  coupling.  Also 
the  equations  are  coupled  in  the  damping  term,  because  the 
system's  two  natural  modes  have  translation  and  rotation 
involved  in  them.  This  model  uses  5  percent  of  the  critical 
damping  coefficient  as  system  damping.  To  evaluate  damping,  a 
modal  analysis  for  the  system  is  performed.  This  analysis  is 
performed  in  Appendix  3.  Substituting  damping  coefficients 
into  Equations  4.3.9a  and  4.3.9b,  the  equations  of  motion  can 
be  expressed  in  matrix  notation  as  Equation  4.3.10: 

CM3  (Yl)  +  [C]  (iLi)  +  [K]  (Yl)  =  {E ' (t) }  (4.3.10) 
where 

{ y ' }  =  Response  Vector  = 


[M] 

=  Mass  Matrix  = 

-mn 

m12" 

-m21 

m22  - 

rci 

=  Damping  Matrix  = 

"c 

11 

C12~ 

c 

21 

C22. 

[K] 

=  Stiffness  Matrix  : 

'kll 

o 

_o 

k22 

{ E '  ( t ) }  =  Seismic  Forcing  Vector  = 


To  simplify  the  explanation,  the  matrices  will  be  redefined  as 
[A]  {y}  +  [B]  (y)  =  (E(t) )  (4.3.11) 

where 


'{e ' (t)  }J 


Equations  4.3.6,  4.3.8,  and  4.3.11  will  be  solved  in  Sections 
6  and  7 . 

4.4  Three  Degree  of  Freedom:  Rotational.  Horizontal, 
and  Vertical  Translation  Response 

4 . 4 . 1  Approach 

This  section  develops  a  three  degree  of  freedom 
mathematical  model  of  a  drydocked  submarine.  The  idealized 
model  of  the  vessel-ctrydock  system  is  shown  in  Figure  4.8.  In 
this  system,  vessel  motion  in  restricted  to  rotation  about  the 
keel,  and  in  horizontal  transverse  and  vertical  translations 
relative  to  the  keel  and  side  pier  supports.  The  approach  for 
determining  vessel-drydock  failure  is  identical  to  the  two 
degree  of  freedom  case. 

The  three  degree  of  freedom  mathematical  model  predicts 
and  analyzes  submarine  rotational,  horizontal  and  vertical 


transitional  responses  caused  by  seismic  ground  accelerations 
normal  to  the  vessel's  longitudinal  axis.  The  rotational 
response  is  in  radians,  the  translational  response  in  inches, 
and  all  three  responses  are  expressed  as  a  function  of  time. 
The  keel  and  side  pier  forces  generated  by  the  responses  are 
compared  to  applicable  pier  failure  criteria  found  in  section 
3  to  check  blocking  stability. 

4.4.2  Three  Degree  of  Freedom  Equation  of  Motion  Model 
As  in  the  two  degree  of  freedom  case,  the  keel  is  the 
origin  for  the  three  degree  of  freedom  system  with  absolute 
rotation,  6,  and  translation  coordinate,  x,  is  the  horizontal 
vessel  motion  relative  to  the  drydock.  The  linear  relative 
translation  coordinate,  y,  is  defined  with  respect  to  vessel 
and  drydock  translation  relative  to  the  ground.  This 
coordinate  y  describes  the  vertical  motion  of  the  vessel 
hull  relative  to  the  drydock  and  the  drydock  blocks.  Defining 
v  as  the  position  of  the  submarine  keel  relative  to  the 
ground,  then  the  following  relationships  hold; 

Y  =  v  -  yg 

v  =  y  +  yg 

and  v=y+yg  (4.4.1) 

Likewise,  u  =  x  +  xg  (4.4.2) 

The  simplified  vessel-block  system  showing  coordinates  (x,  y 
and  ?)  at  rest  and  when  excited  are  shown  in  Figure  4.8  and 
4.9  respectively. 

As  developed  in  Section  4.3.2,  the  first  and  third 
equations  of  motion  in  the  three  degree  of  freedom  model  are 
MX  +  MKG  +  Cxx  +  (2Ksh  +  Kj^JX  =  -Mxg 


(4.4.3) 


Ik6  +  MKG  X  -  MKSy6  +  Ce 

_  _  ..  _  .. 

+  [(B2/2)Ksv  -  WKG]6  =  -MKG  xg  +  MKG  yg6  .  (4.4.4) 

A  balance  of  forces  in  the  y  direction,  as  shown  in 
Figure  4.9,  yields  the  second  equation  of  motion 

1  Fy  =  Mv  =  (Flsv  -  F1(jv)  +  (f2sv  “  F2dv) 

+  ( F  3  Sy  -  F  3  )—  W.  (4.4.5) 


The  above  equations  can  be  simplified  by  the  following: 

1)  The  vessel  weight,  W,  must  be  equal  to  the  summation 
of  static  block  forces  so  that  static  equilibrium  holds  true, 
i.e. 

Flsv  +  f2sv  "*■  f3sv  =  W  (4.4.6) 

2)  The  dynamic  block  forces,  F^y,  as  seen  from  Figure 
4.9  equals  to  a  force  due  to  the  modelled  springs  plus  a 
dissipative  force  due  to  system  damping.  These  dynamic  block 
forces  are  only  a  function  of  vertical  displacement,  y,  and 
velocity,  y.  The  formulation  of  these  forces  is  identical  to 
the  procedure  used  in  Section  4.3.2  and  will  not  be  repeated. 
The  final  form  of  the  dynamic  block  forces,  F^y,  is 

Fdv  ~  (2KSV  +  Kjcy)y  +  Cy  y  (4.4.7) 

where  Ksv  =  side  pier  vertical  stiffness 
Kkv  =  keel  pier  vertical  stiffness 
Cy  =  system  vertical  damping  coefficient. 

Substituting  Equations  4.4.1,  4.4.2,  4.4.6  and  4.4.7  into 
Equation  4.4.5  and  rearranging  components  yields 

My  +  Cyy  +  (2KSV  +  Kkv)y  =  -Myg.  (4.4.8) 

The  above  equation  along  with  Equations  4.4.3  and  4.4.4 


formulate  the  three  equations  of  motion  for  the  three  degree 
of  freedom  model. 

The  non-linear  terms,  MKGyQ  and  MKGyge,  found  in 
Equation  4.4.4  will  be  removed  at  this  time  leaving  a 
linearized  equation  since  8  is  assumed  to  be  small.  The 
effects  of  the  non-linear  terms  will  be  evaluated  in  Section  6 
and  7. 

Now  the  system  of  equations  for  the  three  degree  of 
freedom  is  as  follows; 

mu  x  +  m3_3  8  +  C]_  x  +  x  =  -mu  xg  (4.4.9a) 

*22  y  +  c2  y  +  k22  y  =  -® 22  yg  (4.4.9b) 

ln3  3  6  +  ®3 1  X  +  C3  6  +  IC338  =  -m3^  Xg  (4.4.9c) 

where  m^  =  m22  =  M 

mi 3  =  m3  2  =  MKG 
m33  =  Jk 

ci  =  system  horizontal  damping  coefficient 
c2  =  system  vertical  damping  coefficient 
c3  =  system  rotational  damping  coefficient 

kll  =  2Ksh  +  Kkh 
k22  =  2Ksv  +  kkv 
k33  =  B2/2  Ksv  -  WKG. 

The  first  and  third  equations  of  motions.  Equations  4.4.11a 
and  4.4.11c,  respectively,  are  coupled  in  the  mass  times 
acceleration  term,  which  is  known  as  inertial  coupling.  Also, 
these  equations  are  coupled  in  the  damping  term  because  two  of 
the  system's  three  natural  modes  have  horizontal  translation 
and  rotation  involved  in  them.  The  second  equation  of  motion, 
Equation  4.4.11b,  is  uncoupled  from  the  rest  of  equations. 


Hence,  the  last  natural  mode  of  the  system  only  depends  on 
vertical  translation. 


In  order  to  introduce  damping,  this  model  uses  5  percent 
of  the  critical  damping  coefficient  as  system  damping.  To 
evaluate  damping,  a  modal  analysis  for  the  system  in  question 
must  be  performed.  This  analysis  is  performed  in  Appendix  3  . 
System  damping  coefficients  can  be  found  using  this  modal 
anlaysis  technique,  but  for  now,  the  coefficients  will  be  left 
in  terms  of  c^,  C13,  C22  /  c3i  and  c33*  Substituting  these 
coefficients  into  Equations  4.4.11a,  4.4.11b,  and  4.4.11c,  the 
equations  of  motion  can  be  expressed  in  matrix  notation  as 
Equation  4.4.12: 

CM]  +  [C]  <il)  +  [K]  {£_)  =  (EJL(t)  }  (4.4.10) 


where 

(y f } 


=  Response  Vector  = 


B 


mn 

0 

m 

[M] 

=  Mass  Matrix  = 

0 

m22 

0 

m31 

0 

m 

13 


33 


[K] 


=  Stiffness  Matrix  = 


kn  0 


•22 


o 

o 


{ E ' (t )  }  =  Seismic  Forcing  Vector  =  x 


0 

k33_ 

_mll 

X 

g 

-m22 

_m31 

X 

g 

[C] 


=  Damping  Matrix  = 


c 

ix  = 


C11  0 


'13 


FIGURE  4.2 

Vessel-Drydock  System  Force  Diagram  with  Applied  Quasi-static 

Force 
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Table  5.2  lists  the  vessel  radii,  weight,  distance  from  center 
of  gravity  to  keel,  KG,  moments  of  inertia  about  vessel  center 


of  gravity,  Ic.g.  and  keel,  1^,  and  side  block  transverse 
spacing,  B. 


TABLE  5.2 

Key  Vessel  Parameters 

Hull 


SYSTEM 

Radius 

(inches) 

Weight 

(Kips) 

KG 

(inches) 

^CCT  _ 

(K-in-sec2) 

Ik  , 

(K-m-sec^) 

B 

( inches 

1-5 

198 

16,396 

193 

831,257 

2,411,000 

144 

to 

i 

VD 

252 

37,656 

223 

3,097,535 

7,949,000 

180 

9 

198 

13,624 

193 

691,852 

2,007,000 

138 

10-11 

190 

9,529 

174 

445,592 

1,193,000 

138 

5_. 

2  Block 

Parameters11 

The  models  developed  in  Section  4  require  the  evaluation 
of  the  block's  vertical  and  horizontal  spring  stiffness 
constants,  Kv  and  The  spring  stiffness  is  the  amount  of 

force  required  to  produce  one  unit  of  displacement  in  the 
respective  direction.  In  order  to  facilitate  the  use  of 
Hooke's  law  in  evaluating  spring  constants,  the  assumption  of 
homogeneous  isotropic  behavior  will  be  used  for  the  block 
materials:  wood  and  concrete.  Both  the  side  and  keel  piers 
must  be  analyzed  for  Kv  and  K^. 

5.2.1  Vertical  Stiffness  Kv11 
Using  Hooke's  Law,  the  linear  relationship  between 
vertical  applied  force  and  deformation  can  be  written  as 


%  S 


where  F  =  vertical  force 


=  area  under  applied  force 
=  modulus  of  elasticity 
=  change  in  height 
=  original  height. 


Hence,  the  spring  constant  for  a  given  material  is 

k  =  l  =  ea 

v  h  H  • 


(5.1) 


For  the  eleven  drydock  configurations,  two  types  of 
standard  blocks  (composite  and  timber)  are  used.  A  standard 
composite  block  which  is  composed  of  a  softwood  cap,  a 
hardwood  middle  portion,  and  a  concrete  bottom  can  be  modelled 
as  a  series  of  three  springs  (Figure  5.1).  The  resultant 


block  spring  stiffness  is 


where 


^v  block  —  [  ( V^cap)  +  (i/^oak)  +  (i/^con)]  (5.2) 

where 

kcap  =  softwood  cap  spring  constant 
koak  =  hardwood  middle  spring  constant 
kCOn  =  concrete  bottom  spring  constant. 

All  three  constants  listed  above  are  calculated  from 
Equation  5.1.  A  standard  timber  block  is  composed  of  a  softwood 
cap  and  a  hardwood  body  (Figure  5.2).  Its  block  spring 


constant  is 


^v  block  —  [  ( V^cap)  +  (i/^oak)^  ^  •  (5.3) 

The  total  vertical  stiffness  for  a  drydock  configuration 
can  be  computed  by  multiplying  kv  block  times  the  number  of 
blocks.  The  keel  pier  has  a  total  vetical  stiffness  of 
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P(E2I_E3I3)H3 
+  3E2E3II3 

where  =  height  of  concrete  (27  inches) 

H2  =  height  of  hardwood  (29  inches) 

H3  =  height  of  softwood  cap  (4  inches) 

E^  =  modulus  of  elasticity  of  concrete 
E2  =  modulus  of  elasticity  of  hardwood 
E3  =  modulus  of  elasticity  of  softwood  cap 
I  =  moment  of  inertia  of  block's  cross-section 


I3  =  moment  of  inertia  of  cap's  cross-section. 

For  the  bending  displacement  of  *  timber  block,  Equation  5.4 
still  holds  true  with  the  alterations: 


H;l  =  height  of  hardwood  (27  inches) 

=  modulus  of  elasticity  of  hardwood. 

In  shear,  deformation  can  be  determined  by  modelling  the 
composite  block  as  an  element  subjected  to  shear  stress  at  the 

top  (Figure  5.4).  The  shear  displacement  becomes 

(1+v  )PH,  2  ( 1  +  v  , )  PH_ 

j  _  con  1  ^  wood  2 

s  AE  AE2 


(5.5) 


where  A  =  area  of  cap's  cross  section 

vwood  =  Poisson's  ratio  for  wood  (0.30  is  used) 
vcon  =  Poisson's  ratio  for  concrete  (0.15  is  used). 


By  redefining  H]_  and  Elf  Equation  5.5  yields  the  shear 
displacement  of  a  timber  block. 


The  horizontal  spring  constant  of  any  one  block  can  be 
determined  through  the  Hooke's  Law  relation 

kh  =  P/(dB  +  ds)  (5.6) 

The  total  horizontal  stiffnesses  for  the  keel  and  side  piers 
are 

Kkh  =  kh  *  nk 
and  Ksh  =  kh  *  ns  . 

The  values  for  the  total  area  of  the  rows  of  blocks  for 
each  of  the  eleven  systems  are  listed  in  Table  5.3.  Also  listed 
are  the  corresponding  total  area  of  the  cap  blocks  for  a  row  of 
side  or  keel  blocks.  The  values  of  total  horizontal  and 
vertical  stiffnesses  of  the  keel  and  side  piers  for  the  eleven 
drydock  configurations  are  listed  in  Table  5.4. 


TABLE  5.4 


Total 

Keel  and  Side 

Pier  Stiffness 

in  Pounds 

System 

Kjcv 

Ksv 

Kjch 

1 

50089159 

7545899. 

14018741 

2 

50089159 

3903051. 

14018741 

3 

28493109 

4292468. 

5161684. 

4 

28493109 

2220242 . 

5161684. 

5 

50089159 

2220242 . 

14018741. 

6 

98356895 

52105736 

24330635 

7 

98356895 

34542004 

24330635 

8 

98356895 

26345596 

24330635 

9 

21758374 

8239775. 

3941650. 

10 

16577809 

6765289. 

3072572. 

11 

16577809 

5464271. 

3072572. 

Per  Inch 
KSh 

2612721 

1351407. 

1282216. 

663215.3 

663215.3 

16104611 

10676090 

8142780. 

2683291. 

2139500. 

1728057. 


6.0  VESSEL  RESPONSE  TO  SEISMIC  EXCITATION 


This  portion  of  the  thesis  discusses  methods  of 
determining  the  seismic  response  of  the  vessel -drydock  system 
mathematical  models  developed  in  Section  4.  The  numerical 
analysis  schemes  to  be  implemented  are  a  modal  analysis  method 
(Section  6.1)  which  looks  at  the  linearized  equations  of 
motion,  and  a  fourth  order  Runge-Kutta  numerical  method 
(Section  6.2)  which  analyzes  the  non-linear  aspect  and  also 
verifies  the  linear  solution.  Section  6.3  incorporates  a 
response  spectrum  analysis  to  solve  the  mathematical  models. 
Finally,  Section  6.4  generates  system  response  for  the  eleven 
vessel-drydock  configurations  using  the  quasi-static  force 
method . 


6 . 1  Vessel  Response  Using  Modal  Analysis 
In  order  to  use  the  modal  analysis  method  described 
herein,  only  the  linearized  equations  of  motion  mathematical 
models  can  be  evaluated.  These  are  Equations  4.2.8,  4.3.11 
and  4.4.11  which  correspond  to  the  one,  two,  and  three  degree 
of  freedom  models  respectively.  Since  all  three  linearized 
equations  are  in  the  form 

[A]{y}  +  [B] (y)  =  (E (t) } , 

the  three  degree  of  freedom  model  will  be  analyzed.  The  other 
two  cases  are  just  a  reduced  version  of  the  third. 

To  solve  the  linearized  three  degree  of  freedom 
equations  of  motion  (4.4.11)  in  matrix  notation  and  obtain 
system  responses  as  a  function  of  time,  the  problem  is 
solved  as  an  eigenvalue  problem  with  a  procedure  described  in 
reference  [13].  As  a  first  step,  the  equations  are  decoupled. 


To  accomplish  this,  first  consider  the  free  vibration  problem 
[A]{£}  +  [B]{y}  =  0  (6.1.1) 

Assuming  Equation  6.1.1  has  a  solution  in  the  form 
(y>  =  (V>  eXt  , 

Equation  6.1.1  can  be  expressed  as 

X[A]U)  +  [B]  { =  0.  (6.1.2) 

where  X  =  complex  eigenvalues 
{ 4}  =  modal  vectors. 

Solving  the  matrix  determinant  problem 

I  A[A]  +  [B] |  =0  (6.1.3) 

yields  the  eigenvalues.  In  the  three  degree  of  freedom  case, 
[A]  and  [B]  are  6x6  matrices  hence  the  solution  of  Equation 
6.1.3  will  yield  6  eigenvalues.  These  eigenvalues  will  be 
comprised  of  three  sets  of  complex  conjugates.  Now  the  modal 
vectors  can  be  found  by  substituting  the  eigenvalues  back  into 
Equation  6.1.2  and  solving.  The  modal  vectors  are  then 
combined  to  form  the  modal  matrix  [  where 

U]  =  [  {  i  (W2  («3  K)4  {<i>5  U>6J- 
The  equations  of  motion  can  then  be  decoupled  by  using  the 


modal  matrix, 

i .  e . 

C  ^]T 

[AH40 

-  [ — A. — ] 

and 

[^]T 

[B]  [ip] 

=  C — B —  ] 

where  [ — ■  A _ .] 

and 

[  '  B  — .] 

are  diagonal  matrices 

Once  the  modal  matrix  is  found,  the  forced  vibration 
problem  is  analyzed.  Assume  the  following 
(¥)  =  [■,']  (z) 

and  (z)  =  [^]_1  (y). 


/v; 


A/./  A 


V  \  \ 


A  _*  -■> 
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where  z  is  a  complex  coordinate  system.  Equation  4.4.1  is  now 
expressed  as 

[A]  m  {z}  +  [B]  [*]  {z}  =  (E  (t) } . 

Using  [^]T,  the  above  equation  can  now  be  in  the  form 

r— a— ]  (z)  +  c  “B i  { z }  =  mT{£}  —  {n (t) } . 

This  will  yield  six  uncoupled  equation  of  motion  in  the  form 
aii  zi (t)  +  bii  Zj_(t)  =  Nj_(t),  i  =  1,6 
or  z*i(t)  -  vj_  Zi(t)  =  (1/aii)  Ni(t)  (6.1.4) 

where  vj_  =  -bii/a^i. 

The  solution  of  Equation  6.1.4  is 


zi(t)  =  —  /  e  vi(t-T)  N.(X)  dT,  i  =  1/6.  (6.1.5) 

aii  o 

Once  Equation  6.1.5  is  solved,  the  (y)  matrix  is  found  by  the 
relationship 


(y>  =  [v]  (z). 


Recalling  from  Section  4.4.2  the  terms  of  the  three 
degree  of  freedom  equations  of  motion  (y)  matrix  is 


{ Y)  - 


Hence,  the  desired  vessel-drydock  system  response  can  be  found 
in  the  fourth,  fifth  and  sixth  terms  of  the  (y)  matrix,  i.e. 


Equations  6.2.1a,  6.2.1b  and  6.2.1c  can  be  made  into  an 
equivalent  first  order  system  by  the  following  substitutions 


x  =  R  =  ei(t,x,R,T) 
y  =  S  =  f^t^S) 
e  =  t  =  gx(t, e,y,R,s,T) 

.  mi3  . 

R  +  T  =  e~  (t ,  x  ,  R,  T) 

mll  2 

S  =  f  2  (t ,  y  ,  S ) 

m 

T  +  R  =  g_  (t,  6,y,R,S,T)  . 

m13  2 


The  Runge-Kutta  formulas  [8]  for  this  system  are 

^11  =  hRn  ^11  =  h(e2 (t,xn,Rn,Tn) } 

fen  =  hSn  *12  =  h{ f2 (t,yn, Sn) } 

fell  =  hTn  l13  =  h{g2(t,en,  yn,Rn,Sn,Tn)} 


k 


21  =  h(Rn  +2'fll)  l21  =  h{e2(t,xn  +  2fcll»Rn  +  j£ll'Tn  +  ^£13 


fe22  =  h(Sn  +2^12)  ^22  ~  M^^'Yn  +  2fe12'Sn  +  |^12)} 

1  i  . z 

fe23  =  h(Tn  +2^13)  ^23  =  h{<32(t'en  +  2fe13'yn  +^12' 


Rn  +  5^11»  sn  +  jt-12'  Tn  +  ^13)  > 

fe31  =  h(Rn  +  2^21)^31  =  h{e2(t,xn  +  2fe21'Rn  +  y^21'Tn  ^23)  > 
1  1  ^ 

fc32  =  h(Sn  +  2^22)^32  =  h(f2(t^n  +  2fc22'sn  +  i-^22)  ) 

1  1  ^ 

fc33  =  h(Tn  +  2^23)^33  =  h(92  6n  +  2k23'Yn  +  |fe22'Rn  +  ^21 


1*  1»  % 
sn  +  2^22 ,Tn  +  2^23 ) > 
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fc41 

“  k(Rn 

+  £3l) 

£41 

= 

h{e2(t,xn  + 

^31'  Rn 

+ 

l31^n  +  l33)) 

k42 

=  h(Sn 

+  ^32 ) 

£42 

= 

h{f2(t/yn  + 

k32>  sn 

+ 

l32)) 

k43 

=  h(Tn 

+  ^33 ) 

£43 

= 

h{g2(t,en  + 

k33>  Yn 

+ 

fe32'Rn  +  ^31/ 

sn  + 

l32>  Tn 

+ 

^33)  > 

and 

lastly, 

xn+l  = 

xn  + 

+ 

2fe21  +  2fc31 

+  fe41) 

yn+i  = 

Yn  + 

h  fci2 

+ 

2fe22  +  2  22 

+  k42) 

en+l  = 

en  + 

fc13 

+ 

2fc23  +  2fc33 

+  fc43) 

Rn+1  = 

Rn  + 

^-(£11 

+ 

2^21  +  2^31 

+  £4l) 

sn+l  = 

sn  + 

|(£12 

+ 

2^22  +  2-C32 

+  ^42 ) 

Tn+l  = 

Tn  + 

^(£13 

+ 

2^23  +  2^33 

+  ^43) 

where  h  is  the  time  step.  Hence,  the  desired  vessel-drydock 
system  response  can  be  found  and  system  failure  analyzed  by 
the  iteration  of  the  Runge-Kutta  scheme. 


6.3  Vessel  Response:  Response  Spectrum  Analysis 
The  degrees  of  freedom  mathematical  models  developed  in 
previous  sections  produce  a  time  history  of  the  seismic 
response  of  a  drydocked  vessel.  These  models  use  numerous 
computer  iterations  to  find  their  time  history  in  order  to 
evaluate  the  system’s  maximum  seismic  response.  Another  way 
to  determine  the  maximum  seismic  response  of  a  drydocked 
vessel  is  the  response  spectrum  analysis  method.  Simplicity 
and  fast  results  are  characteristic  of  this  spectrum  method. 
However,  only  maximum  response  is  produced  and  no  time  history 
is  generated. 


The  response  spectra  are  graphs  of  the  maximum  seismic 
response  of  single  degree  of  freedom  systems  over  a  range  of 


system  natural  frequencies.  These  graphs  are  generated  for 
particular  earthquake  and  can  contain  acceleration,  velocity, 
displacement,  or  combination  of  all  three  responses  for  the 
given  system.  Simplified  response  spectra  are  plotted  for 
earthquake  design  analysis.  The  El  Centro  earthquake  response 
spectra  graph,  shown  in  Figure  6.1,  will  be  used  to  confirm 
the  results  of  the  mathematical  models  in  this  thesis.  The 
simplified  spectrum  portion  of  this  graph  is  described  in 
reference  [2]  and  updated  in  reference  [7]. 

The  simplified  El  Centro  earthquake  response  spectra  used 
described  the  maximum  relative  displacement  of  the  system 
shown  in  the  insert  of  Figure  6.1.  The  linear  equation  of 
motion  for  this  simple  system  is 

mu  +  cu  +  ku  =  -m  ysa  fa(t) 

where  u(t)  =  relative  displacement  of  the  mass  with 

respect  to  the  support 

and  ysa  =  maximum  support  acceleration. 

The  natural  frequency  of  the  system  is  given  by 
u.'2  =  k/m 
f  =  u /2~  . 

Using  the  response  spectra  graph  and  the  particular  frequency, 
f,  of  the  system  of  interest,  the  maximum  value  of  the 
relative  displacement,  um,  is  determined.  In  equation  form, 
the  relative  displacements  of  a  single  degree  of  freedom 
system  with  5%  critical  damping  subjected  to  the  El  Centro 
earthquake  ground  motions  are: 
for  to  <  2.24  rad/sec. 


umax  —  1*4  (l^sa) max  —  11*62  inches,  (6.3.1) 
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for  2.24  <  u'  <  12.74, 

umax  =  (Ysa) max/ w  =  26.03/lo  inches,  (6.3.2) 

for  w  >  12.74, 

umax  =  2.6  (ysa)max/w2  =  331. 53/^2  inches.  (6.3.3) 
These  formulas  generated  the  simplified  response  spectrum 
shown  in  Figure  6.2. 

In  order  to  use  Figure  6.2  to  verify  the  one,  two  and 
three  degree  of  freedom  vessel-drydock  models,  some 
manipulation  must  be  performed.  This  is  because  the  linear 
equations  of  motion  of  the  models  do  not  exactly  match  those 
of  the  system  pictured  in  Figure  6.1.  Examining  the  one 
degree  of  freedom  linear  equation  of  motion  (Equation  4.2.7); 
IfcP  +  25lkun6  +[(B2/2)KSV  -  WKG] 6  =  -M  KG  Xg 


e  +  2£w  n6  +  ([(B2/2)Ksv  -  W  KG]/Ik)6  =  - (M  KG/Ik)  xg. 

The  natural  frequency  for  the  above  equation  is  defined  as 
^n2  =  ((B2/2)Ksv  -  W  KG)/Ik. 

Using  Figure  6.2,  this  frequency  gives  a  corresponding  umax. 
Hence,  the  maximum  rotational  response  of  the  system,  ,  is 
then  determined  by  the  relation 

emax  =  M  KG/Ik  umax. 

Table  6.1  lists  the  maximum  rotational  response  for  the  eleven 
configurations  found  in  Section  5. 


»  ,  *  .  *  ,  ■  ,  •  ►  *»  *» 
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TABLE  6.1 


One  Degree  of  Freedom  Vessel-Drydock  System 
Maximum  Response  Using  Response  Spectrum  Analysis 


System 

U3n 

(rad/sec} 

emax 
(rads. ) 

1 

5.580 

0.01594 

2 

3.933 

0.02619 

3 

4.141 

0.02149 

4 

2.870 

0.03100 

5 

2.870 

0.03100 

6 

10.254 

0.00698 

7 

8.327 

0.00860 

8 

7.255 

0.00987 

9 

6.147 

0.01445 

10 

7.253 

0.01299 

11 

6.498 

0.01450 

The  response  spectrum  analysis  to  the  two  and  three 
degree  of  freedom  systems  is  applied  by  means  of  the  modal 
method  described  in  reference  [2].  Concisely,  the  modal 
method  requires  that  the  system  natural  modes  be  determined. 
Then,  the  response  of  the  system  to  a  known  forcing  input  can 
be  developed  by  mode  superposition.  Treating  each  natural 
mode  as  a  single  degree  of  freedom  system  with  its  own  natural 
frequency,  this  method  is  good  for  simplifying  seismic 
response  analysis.  Since  each  natural  mode  acts  as  a  separate 
system,  the  maximum  response  may  be  determined.  Then,  a 
conservative  estimate  for  the  maximum  response  of  the  original 
system  can  be  made  by  adding  up  the  maximum  response  of  each 
mode. 

The  method  of  adding  maximum  modal  responses  gives  an 
upper  bound  for  the  maximum  system  response.  Since  there  is  a 
known  response  spectrum,  the  amount  in  which  each  natural  mode 
contributes  to  the  maximum  response  for  a  given  input  can  be 


determined.  This  is  known  as  the  participation  factor  of  each 
given  mode.  A  full  description  and  derivation  of  partici¬ 
pation  factors  including  implementation  techniques  can  be 
found  in  reference  [2],  and  will  not  be  presented  here.  The 
participation  factor  for  each  natural  mode  in  the  two  and 
three  degree  of  freedom  case  is  calculated  to  confirm  computer 
generated  results  of  maximum  responses  in  Appendix  3 . 

6.4  Vessel  Response;  Quasi-Static  Force  Analysis 


4.1.4,  i.e. 


6  =  Ms/L2  Ksv 


Mg  =  448  *  A  *  KG 


where 

A  =  displacement  in  tons 
L  =  B/2 . 

The  results  of  this  method  are  listed  in  Table  6.2. 


force 

This  analysis 

ill 

hence,  the  quasi- 
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TABLE  6.2 


System  Response  using  Quasi-Static  Force  Analysis 


SYSTEM 


e.  IN  RADIANS 


S  V  s 


7 . 0  COMPUTER  SOLUTIONS 


In  this  section  of  the  thesis,  the  computer  programs 
needed  to  implement  the  mathematical  models  and  the  associated 
seismic  response  are  developed.  Due  to  the  numerous 
iterations  that  is  required  for  the  nodal  analysis  and 
Runge-Kutta  methods  to  effectively  determine  vessel-drydock 
response,  all  computer  work  is  performed  at  the  MIT  Joint 
Computer  Facility  (JCF) .  The  required  computer  programs  are 
written  in  Fortran  77. 

In  order  to  check  the  computer  generated  results,  the 
vessel-drydock  mathematical  models  are  subjected  to  a 
sinusoidal  earthquake  input.  The  advantage  being  that  closed 
form  solutions  of  responses  to  sinusoidal  inputs  can  be 
calculated.  Thus,  closed  form  response  can  be  compared  to  the 
computer  generated  responses. 

Upon  obtaining  the  correct  results  for  a  sinusoidal 
earthquake,  the  El  Centro  earthquake's  acceleration  time 
history  is  applied  to  the  vessel-drydock  system.  The  maximum 
value  of  these  computer  generated  results  is  compared  to  the 
maximum  response  for  the  system  predicted  by  the  response 
spectrum  method  found  in  Section  6.3. 

7 . 1  Modal  Analysis  Solution 

7.1.1  Computer  Program  Development 

As  described  in  Section  6.3,  the  program  approaches  the 
problem  by  first  assembling  matrices  A  and  B.  The  program 
then  uses  the  International  Mathematical  and  Statistical 
Library  (IMSL)  subroutine  EIGZF  resident  on  JCF  to  perform  the 
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system  free  vibration  eigenvalue  problem  on  matrices  A  and  B. 
The  output  of  subroutine  EIGZF  is  the  modal  matrix  [>;■  ]  of  the 
vessel-drydock  system.  Reference  [4]  contains  further 
information  about  subroutine  EIGZF. 

The  equations  of  motion  can  now  be  de-coupled  by  using 
the  modal  matrix.  De-coupling  the  equations  of  motion  leaves 
equations  of  the  form 

an  Zj_(t)  +  bj_j_  Zj_(t)  =  Nj_(t)  ,  i  =  1,  m 
where  m  is  equal  to  twice  the  number  of  the  particular  degree 
of  freedom  model,  i.e.  for  the  one  degree  of  freedom  -  m  =  2. 
The  solution  to  the  above  equation  is 


2i(t)  =  rrr 


-*/ 


-bii/aii (t-T) 


Nj_ ( t )  dx,  i  =  l,m, 


where  a^  and  b^^  represent  the  diagonal  elements  of  the 
diagonalized  A  and  B  matrices,  respectively.  The  integral 


-bii/aii  (t_T) 


Ni(x)  dt 


can  be  evaluated  using  Simpson's  rule,  which  is 

ZAi(t)  =  [(Ai(t  -  2  Ax)  +  Ni(t  -  2  At) ) e- (bii/aii) 2 At 


+  4Ni(t  -  AT)e"(bii/aii)Ar+  Ni(t)] 
where  At  =  data  time  increments 

and  A (t  -  2  At)  =  IA(t)  evaluated  two  time  increments 

previously. 

The  value  Zi(t)  is  evaluated  as 

z  i  ( t )  =  (1/aii)  (  At/3)  (ZAi(t) )  ,  i  =  l,m. 


Wv 


>  /  /  /  i 


To  obtain  the  values  of  the  system  response  for  a 


particular  time,  the  relationship 

{y}  =  [i|>]  {z} 

where  {y}  =  system  response  vector  is  used,  using  the  column 
vector  [z]  for  the  time  at  which  results  are  desired.  The 
values  of  block  forces  caused  by  the  system  response  are 
calculated  every  0.02  seconds  for  the  first  30  seconds  of  the 
earthquake.  Only  the  first  30  seconds  are  needed  due  to  the 
nature  of  the  El  Centro  earthquake  the  seismic  acceleration 
after  30  seconds  are  small  when  compared  to  the  first  30 
seconds.  The  resultant  block  forces  and  vessel  motions  are 
then  checked  against  failure  criteria  established  in  Section 
3.  The  amplitude  of  the  seismic  acceleration  is  reduced  until 
no  block  failure  occurs.  Hence,  the  maximum  earthquake 
acceleration,  in  g's,  that  the  vessel-drvdock  system  can 
withstand  without  failure  can  be  found. 

The  three  degree  of  freedom  vessel  response  computer 
program  is  listed  in  Appendix  1.  Since  all  three  linearized 
models  found  in  Equations  4.2.8,  4.3.11  and  4.4.11  are  in  the 
form 

[A] {y}  +  [B] { y }  =  { E ( t) } , 

the  three  degree  of  freedom  computer  program  is  modified  to 
give  the  one  and  two  degree  model  results. 

7.1.2  Computer  Program  Input 
The  computer  program  listed  in  Appendix  1  evaluates  the 
horizontal  and  vertical  translations  and  rotation  of  the  three 
degree  of  freedom  vessel-drydock  system  in  increments  of  0.02 
seconds.  The  program  requires  the  following  input:  vessel 


weight,  WEIGHT,  in  kips,  keel  to  center  of  gravity  height 
(KG) ,  H,  in  inches,  vessel  mass  moment  of  inertia  about  the 
keel,  IK,  in  kip-in-sec2,  time  increments  of  ground 
acceleration  data,  DTAU,  in  seconds,  side  and  keel  piers 
vertical  and  horizontal  stiffnesses,  Ksv,  Kj^y,  Kg^,  and  Kj^, 
in  kips/inch,  gravitational  constant,  GRAVITY,  in  inch/sec2, 
side  and  keel  pier  dimensions  of  width  and  height,  BASE  SIDE, 
BASE  KEEL,  HTSIDE  and  HTKEEL,  in  inches,  the  coefficient  of 
friction  between  the  block-block  and  ship-block  interfaces,  U1 
and  U2,  horizontal  distance  between  the  center  of  vessel 
contact  with  port  and  starboard  side  piers,  BR,  in  inches, 
vertical/horizontal  ground  acceleration  ratio,  AMP, 
proportional  limits  of  side  and  keel  block  materials,  PLSIDE 
and  PLKEEL,  in  kips/in2,  total  cross-sectional  area  of  side 
and  keel  piers,  SIDE  AREA  and  KEEL  AREA,  in  in  2,  percent 
critical  damping,  ZETA,  vessel  hull  number,  HULL,  and  finally 
the  vessel-drydock  configuration  number,  NSYS.  The  following 
is  where  these  inputs  can  be  found: 

Table  5.1  NSYS  and  HULL 

Table  5.2  WEIGHT,  H,  Ik,  and  BR 

Table  5.3  BASESIDE,  BASEKEEL,  HTSIDE,  SIDEAREA 

and  KEELAREA 

Table  5.4  ^kv'  Ksv,  ^kh  ^sh* 

The  rest  of  the  inputs  are  assumed  to  be 
GRAVITY  =  384  in/sec2 
U1  =  U2  =  0.5 

PLSIDE  =  PLKEEL  =0.7  kips/in2 
and  AMP  =  0.5. 
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The  value  of  0.5  for  the  vertical/horizontal  ground  acceler¬ 
ation  ratio  (AMP)  parameter  is  used  because  no  vertical 
component  of  the  El  Centro  earthquake  acceleration  was 
available  at  the  time  of  this  thesis. 

Damping  coefficients  equal  to  5%  of  critical  damping 
(i.e.  ZETA  =  0.05)  are  assumed  for  the  entire  system.  This 
analysis  for  damping  coefficients  can  be  found  in  Appendix  3 
and  is  implemented  directly  into  the  computer  program. 

Certain  idealized  system  failure  criteria  from  Section  3 
need  to  be  modified  to  accurately  represent  the  vessel-drydock 
system  as  it  actually  exists.  As  stated  in  Section  3.1,  block 
crushing  occurs  when  the  stress  on  the  block  exceeds  the 
block's  proportional  limit.  The  vessel  rests  upon  soft  wood 
caps  placed  on  top  of  the  keel  and  side  piers  when  in  the 
drydock  cradle.  Crushing  of  the  cap  is  not  considered  a 
failure  due  to  its  small  size.  For  the  blocking  system 
considered,  the  generated  stresses  are  transferred  through  the 
cap  to  the  top  drydock  block. 

Also,  a  side  pier  cannot  tip  inboard  because  the  vessel 
hull  is  not  physically  attached  to  the  pier  and  cannot  pull 
the  pier  beyond  it's  upright  position.  The  softwood  cap  on 
top  of  the  side  pier  causes  the  vessel  force  to  be  applied  12 
inches  from  the  inboard  edge.  As  shown  in  Figure  7.3,  this 
causes  the  static  force  vector  to  fall  outside  of  the  middle 
one-third  of  the  pier  base,  as  discussed  in  Section  3.  Since 
the  pier  will  only  tip  in  the  outboard  direction,  the 
resultant  vessel  force  must  fall  in  the  outboard  one-third  of 


Sin  03, T  di. 
a 


and  B (t)  =  /  E (t)  — - 

J  ^t 

o  e 

The  results  of  a  computer  run  with  vessel-drydock  system 
#1  and  sinusoidal  excitation  (Equation  7.1.1)  indicate  that 
the  maximum  rotation  angle  e  is  0.038243  radians  at  a  time  of 
2.02  seconds.  To  confirm  this  value,  Equation  7.1.2  has  to  be 
evaluated 


e(2.02)-ps  [A  (2 . 02 )  sin  o)^  (2.02)  -B(2.02)  coso3d  (2 . 02 )  ] 


Iv<+, 


k^d 


2.02 


-f 


where  A(2.02)  =  /  323.95  sin(3.92x) 


Cu3nT 


2 . 02£to 

n 


cosu)d  t  di 


2.02 


and  b (2 . 02)  = 


/ 


£umx 

323.95  sin  ( 3 . 92t  ) 


sin^d  t  dx 


From  Section  6.3, 

on  =  [((B2/2)KSV  -  W  KG)/IK)V2 
=  5.580  rad/sec 
and  od  =  un(l  -  E2)  x/2 

=  5.573  rad/sec. 


With  these  values  of  natural  and  damped  natural  frequencies, 
the  value  of  A (2. 02)  can  be  evaluated  using  the  relationship 


/ 


eat  sinbt  cosct  dt  =  eat  [-a  sin(b-c) t- (b-c) cos (b-cjt 

2  [a  +(b -c)j 


+  e 


at  [a  sin (b+c) t- (b+c) cos (b+c) t] 


2 [a^+ (b+c) ^ ] 


Substituting  the  values 


t  =  2.02 
a  =  u)  n  =  .279 
b  =  3.92 

and  c  =  =  5.573 

The  value  of  A (2. 02)  becomes 

323.95(e-5636)-1  *  -  .8307  =  -153.1629 
Similarly,  for  B(2.02), 

t 

at  [  (b-c)sin(b-c) t+acos (b-c) t] 

2  [a2+  (b-c) 2 ] 

t 

at  I  (b+c) sin (b+c) t+acos (b+c) t] 
2[a2+(b+c)2] 

0 

with  the  same  values  of  a,  b  and  c  as  before. 

Evaluating,  B(2.02)  becomes 

323.95  (e-5636)-1  *  -.2191  =  -40.3973 

Now, 

1  426979*1970 

0(2.02)  =  [-153>1629*sin(2. 02*5. 573) 

2.54  2411000. *5. 573  -  (-40 . 3973*cos (2 . 02*5 . 573)  ] 

=  .038247  radians. 

Comparing  the  computer  generated  solution  (0.038243 
radians)  with  the  closed  form  solution  (.038247  radians),  it 
can  be  seen  that  the  computer  program  is  calculating  the  one 
degree  of  freedom  rotation  formula  correctly. 

One  other  simple  method  to  check  to  see  if  the  rotation 
displacement  time  history  is  correct  is  to  look  at  the 
waveform.  At  steady  state  condition,  0(t)  will  be  a 
sinusoidal  with  the  same  frequency  as  the  sinusoidal  ground 
acceleration.  An  amplitude  different  and  a  phase  shift  will 


/ 


at 


sinbt  sinct  dt  =  e 


be  present.  The  steady  state  solution  [12]  of  a  equation  of 
motion  in  the  form 


is 


x 


mx  +  cx  +  kx  =  F  sinwt 


_ 1 _ 

V(l-r2)2+(25r)2 


sm 


tan 


-1 


where  r 

oj  n 
r 

S 

In  this  case. 


r 

£ 

|F| 

k 


ayuin 

(k/m)1/2  =  Natural  Frequency 

1/2  c/  v;  km  =  percent  critical  damping. 

5.580  rad/sec 
3.92/5.580  =  .703 
.05 

_  I  *g  I 

M  KG  — =  1051012.9  lb. in 

2 . 54 

75071464  lb  in. 


so  now, 

6  steady  state  =  .027402  sin(3.92t  —  .1381)  radians  (7.1.3) 
Figure  7.1  represents  the  sinusoidal  ground  acceleration 
(Equation  7.1.1)  and  Figure  7.2  which  represents  the  rotation, 
6(t),  of  the  vessel-drydock  system  as  generated  by  the  modal 
analysis  computer  program.  Note  the  same  frequency  between 
Figures  7.1  and  7.2.  The  steady  state  portion  of  Figure  7.2 
is  matched  to  Equation  7.1.3.  In  conclusion  the  model 
analysis  method  correctly  predicts  the  one  degree  of  freedom 
mathematical  model. 


u  TTO » .H.H.W  .W  W.*  H LWAWWJllVMlW.VAllUl^.W’ll*  E»  T»  ^  A»u»'..«  <.“  -_■  ^  .'.■  ^ -»i 


sis 


7.1.4  Computer  Program  Results 

The  modal  analysis  method  can  now  be  used  with  the  El 
Centro  earthquake  data  to  predict  the  system  response.  Eleven 
vessel-drydock  configurations  (Section  5.0)  are  analyzed  for 
maximum  ground  acceleration  that  the  system  can  withstand 
without  block  failure.  Also,  the  first  block  failure  modes 
exhibited  by  the  systems  are  found. 

For  the  eleven  vessel-drydock  systems,  the  linearized 
one,  two  and  three  degree  of  motion  mathematical  model 
responses  are  calculated  and  listed  in  Table  7.3,  7.4,  and 
7.5,  respectively.  These  tables  can  be  found  on  pages 
thru 9-3.  The  program  further  calculates  that  vessel  liftoff 
(Section  3.4)  will  be  the  first  failure  mode  to  occur  when 
ground  accelerations  are  greater  than  the  maximum  permissible 
ground  acceleration  listed  in  Table  7.3,  7.4,  ari  7.5.  The 
maximum  permissible  ground  acceleration  is  based  on  the 
reduction  of  the  El  Centro  acceleration  amplitude  until  no 
failure  occur.  However,  the  waveform  remains  unchanged. 

In  order  to  get  a  visual  picture  of  the  various 
responses,  vessel-drydock  system  #1  is  selected  to  be  plotted. 
Figures  7.4;  7.5a  and  b  and  7.6a,b,  and  c  are  the  appropriate 
responses  for  the  one,  two  and  three  degree  of  freedom  models. 
Only  the  first  thirty  seconds  of  the  earthquake  are  required 
to  give  the  maximum  system  response  as  seen  by  the  plots. 

Note  that  the  relative  horizontal  displacement  (Figure  7.6a) 
and  rotation  (7.6b)  time  histories  of  the  three  degree  of 
freedom  model  are  indentical  to  the  ones  found  in  the  two 
degree  of  freedom  model.  This  is  due  to  the  uncoupled 


I 


vertical  equation  of  motin  found  in  the  linearized  three 
degree  of  freedom  case. 

As  a  check  for  the  computer  generated  results,  the  system 
maximum  responses  due  to  full  magnitude  El  Centro  earthquake 
accelerations  are  compared  to  the  system  maximum  responses 
determined  by  response  spectrum  analysis  in  Table  7.1  and  7.2. 
These  tables  are  generated  by  combining  the  computer  results 
of  various  models  and  Tables  6.1  (Section  6.3)  and  A. 3. 3 
(Appendix  3) .  Since  the  response  spectrum  method  is  an 
approximate  one  (Section  6.3),  the  maximum  20%  error  between 
its  results  and  the  computer  results  in  the  one  degree  of 
freedom  case  is  reasonable.  Thus,  the  one  degree  program  is 
generating  accurate  results. 

In  the  two  degree  of  freedom  case,  the  response  spectrum 
method  predicts  the  worst  case  if  the  system's  two  natural 
modes  reach  a  maximum  at  the  same  time,  and  add  together. 

This  happened  in  the  rotational  case,  but  only  twice  in  the 
relative  displacement  case.  In  the  majority  of  the  systems, 
the  computer  generated  rotational  response  actually  exceeded 
the  rotational  response  predicted  by  modal  analysis.  The 
maximum  35%  difference  between  the  two  results  can  be 
accounted  for  because  of  the  approximate,  simplified  response 
spectrum  mentioned  earlier.  The  three  degree  of  freedom  case 
has  the  same  order  of  error  as  found  in  the  two  degree  case  as 
shown  in  Table  7.2.  Hence,  all  mathematical  models  perform 
satisfactory  using  the  modal  analysis  computer  program. 


TABLE  7.1 


One  Degree  of  Freedom  Equation  of  Motion  Response  Comparison 
Modal  Analysis  Computer  Results  of  Linear  System 


System 


I  6  I 

radians 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


0.01474 

0.01618 

0.01436 

0.02966 

0.02966 

0.00832 

0.0115 

0.01712 

0.01512 

0.01714 


Response  Spectrum  Method 


System 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


Results  of  Linear  System 


radians 

0.01594 
C. 02619 
0.02149 
0.031 
0.031 
0.00698 
0.0068 
0.00987 
0.01445 
0.01299 
0.0145 


TABLE  7.2 


Three  Degree  of  Freedom  Equation  of  Motion  Response  Comparison 
Modal  Analysis  Computer  Results  of  Linear  System 


1  x  1 

1  y  1 

1  6  1 

Svstem 

( inches} 

( inches} 

( radians} 

1 

0.27971 

0.07374 

0.01458 

2 

0.22661 

0.10948 

0.01636 

3 

0.53288 

0.15395 

0.01546 

4 

0.58122 

0.18352 

0.03130 

5 

0.24721 

0.11703 

0.02992 

6 

0.36179 

0.06123 

0.00814 

7 

0.39118 

0.06870 

0.00879 

8 

0.39889 

0.07140 

0.01154 

9 

0.61065 

0.11749 

0.01587 

10 

0.57808 

0.11970 

0.01547 

11 

0.61730 

0.11676 

0.01734 

Response 

Spectrum  Method  Results 

of  Linear  System 

1  x  | 

1  y  1 

1  6  1 

Svstem 

( inches} 

( inches} 

( radians} 

1 

0.33275 

0.10859 

0.01646 

2 

0.27141 

0.12225 

0.02348 

3 

0.65142 

0.19089 

0.02350 

4 

0.53788 

0.21491 

0.02350 

5 

0.21468 

0.12980 

0.03202 

6 

0.45271 

0.08025 

0.00683 

7 

0.46833 

0.09708 

0.00866 

8 

0.46597 

0.10762 

0.01009 

9 

0.64366 

0.15381 

0.01552 

10 

0.65864 

0.13662 

0.01342 

11 

0.67467 

0.14955 

0.01552 

NOTE:  Two  Degree  of  Freedom  System  Maximum  Response 

correspond  to  x  and  9 . 


7 . 2  Fourth  Order  Runqe-Kutta  Analavsis  Solution 
7.2.1  Computer  Program  Development 
As  described  in  Section  6.2,  the  Runge-Kutta  scheme  uses 
a  set  of  first  order  differential  equations  to  describe  the 
second  order  vessel-drydock  system.  A  simple  example  of  this 
procedure  is  done  on  the  following  equation 
mx  +  cx  +  kx  =  f(t)  +  g(t)x 

let  r  =  x  and  rearranging  yields  two  first  order  equations, 
x  =  r 

r  *  -5  r  -  |  x  +  £(t)  +  g(t)x 

Now,  the  first  order  system  can  be  implemented  into  a  computer 
program. 

Briefly,  the  fourth  order  Runge-Kutta  numerical  scheme  is 
one  that  makes  use  of  predicted  velocity  and  acceleration  at  a 
given  time  to  compute  the  displacement  and  velocity  for  the 
next  time  increment.  This  can  be  done  by  the  basic 
relationships 

xn+i  =  xn  +  xn  At 
and  xn+1  =  xn  +  xn  At 

However,  x  and  x  are  not  found  by  differentiation  but 
through  the  manipulations  of  the  second  order  equation  into  a 
set  of  first  order  equations.  This  can  be  done  by  procedure 
shown  in  the  previous  paragraph. 

The  first  step  in  the  computer  program  development  is  to 


place  the  set  of  first  order  differential  equations  into  the 
Runge-Kutta  formulas  shown  in  Section  6.2.  The  evaluation  of 
these  formulas  will  predict  the  velocity  and  acceleration 


that  will  be  used  in  the  next  time  increment.  Now,  the  system 
response,  both  displacements  and  associated  velocities,  can  be 
calculated  using  the  previous  information.  This  process  is 
marched  through  time  steps  in  our  to  produce  a  time  history  of 
the  system  response. 

Now,  the  values  of  block  forces  caused  by  the  system 
response  are  calculated  every  0.01  seconds  for  the  first  30 
seconds  of  the  earthquake.  The  first  30  seconds  are  only 
required  to  analyze  for  maximum  seismic  response.  The 
resultant  block  forces  and  vessel  motions  are  then  checked 
against  failure  criteria  found  in  Section  3.  The  amplitude  of 
the  seismic  acceleration  is  reduced  until  no  block  failure 
occurs.  Like  in  the  modal  analysis  method  (Section  7.1),  the 
maximum  earthquake  acceleration,  in  g's,  that  the 
vessel -drydock  system  can  withstand  without  failure  can  be 
found. 

The  Runge-Kutta  analysis  gives  the  ability  to  analyze 
both  the  linear  and  non-linear  equations  of  motion  for  the 
one,  two  and  three  degree  of  motion  mathematical  models.  The 
non-linear  three  degree  of  freedom  vessel  response  computer 
program  is  listed  in  Appendix  2.  The  program  is  developed 
around  Equations  4.4.5,  4.4.6  and  4.4.8  and  those  found  in 
Section  6.2.  Since  all  linear  and  non-linear  models  developed 
in  this  thesis  can  be  arranged  into  the  Runge-Kutta  formulas, 
the  non-linear  three  degree  of  freedom  comptuer  program  is 
modified  to  give  desired  model  results. 


The  computer  program  listed  in  Appendix  2  evaluates  the 


horizontal  and  vertical  translations  and  rotation  of  the 
non-linear  three  degree  of  freedom  vessel-drydock  system  in 
increments  of  0.01  seconds.  The  program  requires  the 
identical  input  found  in  Section  7.1.2  for  the  modal  analysis 
method.  Hence,  the  program  inputs  needs  not  be  listed  again. 
Like  in  the  previous  case,  Appendix  3  provides  the  damping 
coefficients  to  be  implemented  into  the  computer  program. 

7.2.3  Computer  Program  Testing 

In  order  to  validate  the  fourth  order  Runge-Kutta 
numerical  scheme,  the  linear  one  degree  of  freedom 
mathematical  model  is  checked  with  a  set  of  sinusoidal  ground 
accelerations  for  input.  The  same  input  ground  motion  that 
is  used  to  validate  the  modal  analysis  is  implemented  in  this 
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E  ( " )  =  xg  =  329.95  sin (3 . 92t) .  (7.2.1) 

This  sample  function  is  selected  because  its  magnitude  and 
frequency  closely  match  that  of  the  maximum  acceleration 
portion  of  the  El  Centro  earthquake.  Also,  vessel-drydock 
system  configuration  #1  is  chosen  for  this  initial  analysis. 
The  results  of  the  linearized  one  degree  of  freedom  case  will 
be  compared  to  the  modal  analysis  and  closed  form  solutions. 

The  results  of  a  computer  run  with  vessel-drydock  system 
#1  and  sinusoidal  excitation  (Equation  7.2.1)  indicate  that 
the  maximum  rotation  angle  G  is  0.038236  radians  at  a  time  of 
2.02  seconds.  Comparing  this  the  computer  generated  solution 
with  the  modal  analysis  solution  (0.038243  radians)  and  the 


closed  form  solution  (0.038247  radians)  found  in  Section 
7.1.3,  it  can  be  seen  that  the  computer  program  is  calculating 
the  linearized  one  degree  of  freedom  rotation  formula 
correctly.  Also,  Figure  7.7  which  represents  the  rotation,  6 
(t) ,  of  the  vessel-drydock  system  as  generated  by  the 
Runge-Kutta  computer  program  is  essentially  identical  to 
Figure  7.2  generated  by  the  modal  analysis  computer  program. 

In  conclusion,  the  fourth  order  Runge-Kutta  method 
correctly  predicts  the  linearized  one  degree  of  freedom 
mathematical  model.  Now,  the  computer  program  will  be  used 
for  the  linear  and  non-linear  one,  two  and  three  degree  of 
freedom  systems. 


7.2.4  Computer  Program  Results 
The  fourth  order  Runge-Kutta  analysis  can  be  implemented 
with  the  El  Centro  earthquake  data  to  predict  the  system 
response.  Eleven  vessel-drydock  configurations  (Section  5.0) 
are  analyzed  for  maximum  ground  acceleration  that  the  ystem 
can  withstand  without  block  failure.  Also,  the  first  block 
failure  modes  exhibited  by  the  systems  are  found. 

For  the  eleven  vessel-drydock  systems,  the  linear  and 
non-linear  one,  two  and  three  degree  of  motion  mathematical 
model  responses  are  calculated  and  listed  in  Table  7.3,  7.4 
and  7.5,  respectively.  These  tables  can  be  found  on  pages 
thru  93.  The  program  further  calculates  that  vessel  liftoff 
(Section  3.4)  will  be  the  first  failure  mode  to  occur  when 
ground  accelerations  are  greater  than  the  maximum  permissible 
ground  acceleration  listed  in  Table  7.3,  7.4  and  7.5. 


In  order  to  get  a  visual  picture  of  the  various  linear 
and  non-linear  responses,  vessel-drydock  system  #1  is  selected 
to  be  plotted.  Figures  7.8,  7.9  a  and  b,  7.10  a,  b  and  c  are 
the  appropriate  response  for  the  linear  one,  two  and  three 
degree  of  freedom  models  and  Figures  7.11,  7.12a  and  b  and 
7.13  a,  b  and  c  are  for  the  non-linear  models.  Note  that 
these  figures  are  almost  identical  to  Figures  7.3,  7.4a  and  b 
and  7.5  a,  b  and  c  generated  by  the  modal  analysis  method 
(Section  7.1.5).  This  reinforces  the  equivalency  of  the 
results  found  by  the  two  numerical  shemes.  This  can  also  be 
seen  by  comparing  the  system  maximum  responses  of  the 
Runge-Kutta  method  to  that  of  the  modal  analysis  method.  Both 
are  listed  in  Table  7.3,  7.4  and  7.5.  The  maximum  error  for 
any  given  case  is  5%.  Hence,  all  mathematical  models  perform 
satisfactory  using  the  fourth  order  Runge-Kutta  computer 


program. 


TABLE  7.3 

Modal  Analysis  Method  Results  of  Linear  System 


System 


e  max 
radians^ 

Time 

(sec} 

Maximum  Acce 
(g's) 

0.003390 

6.64 

0.0759 

0.003883 

9.28 

0.0792 

0.006033 

9.22 

0.1368 

0.006822 

8 . 68 

0.0759 

0.004513 

8 . 68 

0.0462 

0.001999 

5.36 

0.0792 

0.002422 

5.84 

0.0924 

0.002760 

8.98 

0.0792 

0.005138 

7.54 

0.0990 

0.004537 

8.98 

0.0990 

0.004972 

7.92 

0.0957 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Linear  System 


Time 

(sec) 


Maximum  Acceleration 

_ JtJls) _ 


TABLE  7.4 

Two  Degree  of  Freedom  Equation  of  Motion  Response 
Modal  Analysis  Method  Results  of  Linear  System 


System 


xmax 

inches)  ( 

Time 

seconds) 

e 

max 

(radians) 

Time 

( seconds) 

Maximum 
Acceleration 
(e' s) 

0.064335 

8.38 

0.003355 

9.02 

0.0759 

0.063873 

5.38 

0.003927 

9.26 

0.0792 

0.207826 

5.42 

0.006033 

9.24 

0.1287 

0.127868 

5.72 

0.006887 

8.66 

0.0726 

0.032138 

5.80 

0.003891 

8.72 

0.0429 

0.090448 

5.78 

0.002037 

5.36 

0.0825 

0.109531 

5.46 

0.002463 

5.86 

0.0924 

0.095734 

5.98 

0.002770 

7.82 

0.0759 

0.195410 

6.06 

0.005080 

7.60 

0.1056 

0.167644 

6.04 

0.004487 

7.84 

0.0957 

0.172847 

6.08 

0.004856 

7.54 

0.0924 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Linear  System 


Maximum 


System 


max 

finches)  ( 

Time 

'seconds ) 

max 

( radians ) 

Time 

( seconds ) 

Accelerat 
(e’ s) 

0.064566 

8.38 

0.003358 

9.02 

0.0759 

0.064934 

5.38 

0.003927 

9.26 

0.0792 

0.207015 

5.56 

0.006034 

9.24 

0.1287 

0.127317 

5.72 

0.006890 

8.66 

0.0726 

0.032442 

5.80 

0.003891 

8.72 

0.0429 

0.090984 

5.78 

0.002048 

5.36 

0.0825 

0.110229 

5.46 

0.002466 

5.86 

0.0924 

0.091996 

5.98 

0.002656 

7.81 

0.0759 

0.194354 

6.06 

0.005074 

7.59 

0.1056 

0.166760 

6.03 

0.004493 

7.83 

0.0957 

0.174626 

6.07 

0.004860 

7.55 

0.0924 

4th  Order  Runge-Kutta  Numerical  Method  Results  of  Non-Linear  System 


System 


^max  Time  ®max 

inches)  (seconds)  (radians 


Maximum 

Time  Acceleration 
seconds) 


1 

0.064597 

8.38 

0.003358 

0.0759 

2 

0.064839 

5.38 

0.003926 

0.0792 

3 

0.206443 

5.56 

0.006035 

9.2 

4 

0.1287 

4 

0.127281 

5.72 

0.006892 

8.6 

6 

0.0726 

5 

0.032442 

5.80 

0.003891 

8.7 

2 

0.0429 

6 

0.090984 

5.78 

0.002048 

5.3 

6 

0.0825 

7 

0.110274 

5.46 

0.002467 

5.8 

6 

0.0924 

8 

0.092031 

5.98 

0.002656 

7.8 

1 

0.0759 

9 

0.194361 

6.06 

0.005075 

7.5 

9 

0.1056 

0 

0.166763 

6.03 

0.004493 

7.8 

3 

0.0957 

1 

0.174636 

6.07 

0.004859 

7.5 

5 

0.0924 

TABLE  7.5 


Three  Degree  of  Freedom  Equation  of  Motion  Response 
Modal  Analysis  Method  Results  of  Linear  System 


xmax 

A 

Maximum 

Time 

y  max 

Time 

max 

Time 

Acceleration 

System 

( inches ) 

( inches ) 

(radians ) 

■HHH 

(z' s) 

1 

0.067132 

8.38 

0.017698 

5.68 

0.003501 

9.02 

0.0792 

2 

0.061211 

5.38 

0.025181 

8.06 

0.003765 

9.26 

0.0759 

3 

0.207826 

5.42 

0.060041 

5.64 

0.006033 

9.24 

0.1287 

4 

0.122056 

5.72 

0.038540 

5.64 

0.006574 

8.66 

0.0693 

5 

0.032138 

5.70 

0.015215 

5.38 

0.003891 

8.72 

0.0429 

6 

0.090448 

5.78 

0.019865 

5.32 

0.002037 

5.36 

0.0825 

7 

0.105619 

5.46 

0.018550 

6.68 

0.002375 

5.86 

0.0891 

8 

0.091745 

5.98 

0.016424 

5.36 

0.002654 

7.82 

0.0759 

9 

0.201516 

6.06 

0.038773 

5.80 

0.005238 

7.60 

0.1089 

10 

0.161863 

6.04 

0.033516 

5.38 

0.004332 

7.84 

0.0924 

11 

0.179020 

6.08 

0.033861 

5.40 

0.005030 

7.54 

0.0957 

4  th 

Order  Runge -Kutta 

Numerical 

Method 

Results  of 

Linear 

System 

X 

max 

6 

Maximum 

Time 

^max 

Time 

max 

Time 

Acceleration 

System 

( inches) 

( sec . ) 

(  inches') 

( sec . ) 

(radians') 

( sec . ) 

(e’  s') 

1 

0.067373 

8.38 

0.017858 

5.37 

0.003504 

9.02 

0.0792 

2 

0.062228 

5.39 

0.025102 

8.06 

0.003761 

9.27 

0.0759 

3 

0.207015 

5.43 

0.059736 

5.64 

0.006304 

9.25 

0.1287 

4 

0.121529 

5.73 

0.038651 

5.65 

0.006577 

8.67 

0.0693 

5 

0.032442 

5.81 

0.014985 

5.39 

0.003891 

8.79 

0.0429 

6 

0.094623 

5.78 

0.019463 

5.77 

0.002130 

5.37 

0.0858 

7 

0.106292 

5.47 

0.018152 

6.69 

0.002378 

5.87 

0.0891 

8 

0.091996 

5.99 

0.016558 

5.37 

0.002656 

7.82 

0.0759 

9 

0.200428 

6.06 

0.039205 

5.40 

0.005233 

7.60 

0.1089 

10 

0.161009 

6.04 

0.033887 

5.39 

0.004338 

7.84 

0.0924 

11 

0.180863 

6.08 

0.034608 

5.40 

0.005033 

7.55 

0.0957 

4th  Order  Runge 

-Kutta  Numerical  Method  Results  of  Non-linear  System 

X 

max 

®  max 

Maximum 

Time 

-^max 

Time 

Time 

Acceleration 

System  (inches)  (sec.)  (inches’)  (sec.)  (radians')  (sec.) _ (z' s) 


1 

0.067283 

8.38 

0.017858 

5.37 

0.003507 

9.02 

0.0792 

2 

0.062050 

5.39 

0.025102 

8.06 

0.003763 

9.27 

0.0759 

3 

0.206406 

5.57 

0.059736 

5.64 

0.006037 

9.25 

0.1287 

4 

0.121484 

5.73 

0.038651 

5.65 

0.006580 

8.67 

0.0693 

5 

0.032443 

5.81 

0.014985 

5.39 

0.003891 

8.79 

0.0429 

6 

0.094523 

5.78 

0.019463 

5.77 

0.002130 

5.37 

0.0858 

7 

0.106348 

5.47 

0.018152 

6.69 

0.002379 

5.87 

0.0891 

8 

0.091968 

5.99 

0.016558 

5.37 

0.002656 

7.82 

0.0759 

9 

0.200131 

6.06 

0.039205 

5.40 

0.005234 

7.61 

0.1089 

10 

0.160892 

6.04 

0.033887 

5.39 

0.004337 

7.84 

0.0924 

11 

0.180674 

6.08 

0.034608 

5.40 

0.005032 

7.55 

0.0957 

Discussion  of  the  Models  and  Associated  Results 


Both  the  modal  analysis  and  the  Runge-Kutta  numerical 
methods  give  satisfactory  results  in  determining  the  system 
maximum  response  in  the  various  models  and  configurations. 

The  maximum  error  in  the  comparison  of  the  two  methods  for 
any  given  run  is  5%.  This  is  excellent  when  considering  that 
each  method  has  completely  different  basic  concepts  behind 
them.  The  modal  anlysis  method  runs  faster  to  obtain  the 
result  due  to  less  computer  operations.  But,  the  Runge-Kutta 
method  is  much  easier  to  program  into  the  computer  and  is  more 
flexible.  In  all,  both  methods  can  be  used  to  predict  the 
seismic  response  of  the  linearized  one,  two  and  three  degree 
of  freedom  models. 

In  the  comparison  of  the  one,  two  and  three  degree  of 
freedom  responses,  certain  conclusions  can  be  drawn.  The  one 
degree  of  freedom  model  need  only  be  used  if  system  maximum 
rotational  response  is  the  main  concern  to  drydock  system 
design.  This  model  gives  results  sufficiently  close  to  the 
rotation  generated  by  the  other  two  models.  Since  vessel 
liftoff  of  the  side  block  is  the  first  block  failure  to  be 
exhibited  by  the  vessel-drydock  system,  the  one  degree  of 
freedom  model  adequately  predicts  the  maximum  earthquake 
amplitude  the  system  can  withstand.  Most  of  the  vertical 
displacement  at  the  side  pier  is  due  to  rotation  and  not  of 
relative  vertical  displacement  as  found  in  the  three  degree 
model.  Hence,  the  one  degree  of  freedom  system  response 
should  be  used  to  determine  vessel  roation  and  maximum 
earthquake  amplitude. 


The  three  degree  of  freedom  system  results  should  be  used 


for  a  parameter  study  on  block  failure  modes  and  for  block 
design.  This  is  because  only  the  three  degree  model  produces 
all  required  displacements  for  estimation  of  generated  forces 
in  both  keel  and  side  blocks.  However,  when  the  relative 
vertical  translations  generated  are  small,  the  two  degree 
model  can  be  used  in  block  design, .  The  one  degree  model  is 
truely  inadequate  for  block  design  since  crushing  and  vessel 
liftoff  of  the  side  blocks  are  the  only  failure  modes. 

The  non-linear  equations  of  motion  mathematical  models 
needs  to  be  addressed.  As  shown  in  Table  7.3,  7.4  and  7.5,  the 
non-linearity  of  th ■».  system  caused  less  than  1%  difference 
between  the  linear  and  non-linear  responses.  Under  the  eleven 
respresentative  vessel-drydock  configurations,  the  non-linear 
terms  (ie.  y6  and  ygG)  should  be  removed  from  the  equations  of 
motion.  This  linearization  of  the  equations  of  motion  will  be 
verified  in  the  proceeding  discussion. 

The  next  step  in  the  studying  of  the  non-linearity  of  the 
system  is  to  determine,  when  the  non-linear  terms  mentioned 
earlier  should  be  included  in  the  equations  of  motion.  The 
non-linear  equation  of  motion  for  the  one  degree  of  freedom 
model  is 


-M  KG  ••  M  KG  ..  „ 

i7~  xg  f  T T  yg  6 


Let  Xg  =  a  sin  t 
YG  =  ' xg 
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M  KG  a/Ik  =  P 
A2  =  ki/Ik 

and  2  5 A  =  Cq /Ik* 

Substituting  and  rearranging  the  Equation  of  motion  gives 

**  • 

0  +  2£A0  +  A26  =  -6  sinut(l  -yQ)  .  (7.3.1) 

where  A  *  system  natural  frequency 

6  =  normalized  maximum  amplitude 
y  =  ratio  of  vertical/horizontal  amplitude. 

Equation  7.3.1  not  only  represents  the  vessel-drydock  system 
but  any  other  system  than  can  be  reduced  to  this  form.  A 
parameter  study  is  performed  using  \  3  and  Y  to  determine 
whether  or  not  to  use  non-linear  terms  in  the  equations  of 
motion.  The  selection  criterion  will  be  when  the  non-linear 
response  (ie.  Y  /  0)  is  different  to  the  linear  response  (ie  Y 
=  0)  by  10%.  If  for  a  given  A ,  8  and  Y  this  criterion  is 
exceeded,  then  the  non-linear  term,Y6,  should  be  included.  If 
not,  exclude  them.  For  this  study,  let  u  «=  3.92  rad/sec  and  5 
=  0,05  to  emulate  the  sinusoidal  excitated  system  with  5% 
damping  discussed  in  Section  7.1.3.  Using  a  Runge-Kutta 
numerical  scheme,  Figures  7.14  a,  b,  c,  d  and  e  are  generated 
by  varying  the  three  parameters.  Note  that  a  negative  Y 
indicates  a  phase  shift  of  180°  between  the  vertical  and 
horizontal  sinusoidal  excitation.  The  linearity  selection 
criterion  of  10%  difference  holds  true  for  the  region  marked 
"N-L”  for  non-linear  for  a  given  6 .  in  the  case  of  the  one 
degree  of  freedom  vessel -drydock  systems,  the  linearizing  of 
the  equations  of  motion  can  be  shown  on  Figure  7.14b  since  3  = 
0 . 5 ,Y *  0.5  and  A  ranges  from  2.870  rad/sec  to  10.254  rad/sec 


(Table  6.1) . 


In  order  to  verify  the  linearity  selection  criterion,  a 
linear  point  selected  off  of  Figure  7.14b  of  B  =  0.5,  Y  =  5.58 
andY  =  20.0  is  implemented  into  the  Runge-Kutta  computer 
program  in  Section  7.2.  The  nonlinear  response  equals 
0.015839  radians  compared  to  the  linear  response  of  0.014741 
radians  when  using  the  El  Centro  earthquake  input.  This  is  a 
difference  of  7.5%  which  is  less  than  10%.  The  linearity 
selection  criterion  holds  true.  Now,  a  non-linear  point  of  B 
=0.5,  X  =  8.327  and  y  =  15.0  is  implemented  into  the 
Runge-Kutta  scheme.  The  non-linear  and  linear  responses  are 
.009652  and  ,008664,  respectively.  This  is  a  difference  of 
11.4%  which  is  greater  than  10%.  The  linearity  selection 
criterion  is  validated.  Thus,  vessel-drydock  system  should  be 
checked  before  selection  of  linear  or  non-linear  models. 

As  mentioned  earlier.  Equation  7.3.1  can  be  used  to 
validate  systems  other  than  the  vessel-drydock  systems 
subjected  to  sinusoidal  ground  acceleration  whose  frequency 
matches  that  of  the  El  Centro  earthquake,  ie.  ^  =  3.92 
rad/sec.  this  can  be  done  by  modifying  the  appropriate 
linearity  selection  criteria  figure  with 
Xj  =  X/3.92 

and  replotting  the  figure.  This  essentially  changes  Equation 
7.3.1  into  a  r on-dimensional  of 

0]_  +  2  C  +  X^2  s =  -  (3^  sinti  (1  —  Y^6]_) 


where 
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FIGURE  7.2 

Vessel  Seismic  Response  in  Rotation  using  Modal  Analysis  with 
Sinusoidal  Ground  Acceleration 
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FIGURE  7.4 

Vessel  Seismic  Response  in  Rotation  for  One 
Model  using  Modal  Analysis 
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FIGURE  7.6b 

Vessel  Seismic  Response  in  Relative  Vertical  Translation  for 
Three  Degree  of  Freedom  Model  using  Modal  Analysis 
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FIGURE  7.7 

Vessel  Seismic  Response  in  Rotation  using  Runge-Kutta 
Numerical  Analysis  with  Sinusoidal  Ground  Acceleration 
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FIGURE  7.9a 

Vessel  Seismic  Response  in  Relative  Horizontal  Translation  for 
Two  Degree  of  Freedom  using  Runge-Kutta  Analysis 


RELATIVE  HORIZONTAL 


0  0  10  15  20  20  30 


TIME  IN  SECONDS 


FIGURE  7.10a 

Vessel  Seismic  Response  in  Relative  Horizontal  Translation  for 
Three  Degree  of  Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.10c 

Vessel  Seismic  Response  in  Rotation  for  Three  Degree  of 
Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.11 

Vessel  Seismic  Response  in  Rotation  for  Non-Linear  One  Degre 
of  Freedom  Model  using  Runge-Kutta  Analysis 
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FIGURE  7.13a 

Vessel  Seismic  Response  in  Relative  Horizontal  Translation 
for  Non-linear  Three  Degree  of  Freedom  Model  using  Runge- 

Kutta  Analysis 
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FIGURE  7.13b 

Vessel  Seismic  Response  in  Relative  Vertical  Translation  for 
Non-linear  Three  Degree  of  Freedom  Model  using  Runge-Kutta 

Analysis 
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FIGURE  7.14c 

Linearity  Selection  Criteria,  5%  Damping,  Beta 
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8.0  SUMMARY.  CONCLUSIONS.  AND  RECOMMENDATIONS  FOR  FURTHER 

STUDY 

8 . 1  Summary 

This  thesis  investigates  four  models  of  the  vessel- 
drydock  system  for  seismic  analysis.  These  models  are 
implemented  to  determine  vessel  response  to  seismic  loading 
for  eleven  representative  submarine -drydock  systems.  Using 
Equation  4.1.4,  the  quasi-static  force  model  generates  various 
vessel  responses  listed  on  Tabe  6.2.  For  example,  the 
submarine  in  system  #1  will  rotate  about  its  keel  0.00809  when 
subjected  to  a  strong  earthquake.  The  one  degree  of  freedom 
model  predicts  a  maximum  vessel  rotation  for  system  #1  of 
0.01471  radians  in  response  to  the  El  Centro  earthquake  ground 
acceleration  history.  This  rotation  of  0.01471  radians  is 
verified  by  the  two  and  three  degree  of  freedom  models.  See 
Tables  7 . 2  and  7.3.  This  maximum  seismic  response  of  system 
#1  can  also  be  confirmed  by  the  response  spectrum  method  as 
described  in  Section  6.3  and  reference  [1],  The  response 
spectrum  method  gives  the  rotation  as  0.01594  radians.  All 
eleven  systems  have  similar  differences  between  the  quasi¬ 
static  force  model  and  the  various  degree  of  freedom  models. 

More  important  than  the  maximum  vessel  response  due  to 
the  El  Centro  earthquake  of  .33g  magnitude  is  the  maximum 
permissible  earthquake  acceleration  that  a  vessel -drydock  can 
withstand  without  failure.  The  maximum  permissible 
accelerations  and  associated  vessel  responses  for  the  one,  two 
and  three  degree  of  freedom  are  listed  in  Table  7.3,  7.4  and 
7.5,  respectively.  For  example,  the  one  degree  of  freedom 


v’-.'V 


model  for  system  #1  predicts  a  maximum  permissible 
acceleration  of  0.0759  g's  with  maximum  vessel  rotation  of 
0.00339  radians.  The  two  degree  of  freedom  determines  a 
maximum  permissible  acceleration  equalling  .0759  g's  with 
maximum  vessel  rotation  and  relative  horizontal  displacement 
of  0.003358  radians  and  0.064566  inches,  respectively. 

Finally,  the  maximum  permissible  acceleration  of  0.0792  g's 
with  associated  vessel  responses  (rotation,  relative 
horizontal  and  vertical  displacements)  of  0.003504  radians, 
0.067373  inches,  and  0.017858  inches.  The  slightly  higher 
maximum  permissible  acceleration  found  in  three  degree 
model  is  due  to  the  beneficial  effect  of  the  vertical 
displacement  delaying  failure  of  the  system.  However,  as 
shown  by  the  other  systems,  this  is  not  always  the  case.  In 
summary,  the  range  of  magnitudes  known  as  maximum  permissible 
accelerations  that  the  eleven  systems  could  withstand  without 
liftoff  is  13%  to  42%  of  the  magnitude  of  the  El  Centro 
earthquake  depending  on  system  and  model  used. 

8 . 2  Conclusions 

The  quasi-static  method  currently  used  by  the  U.S.  Navy 
for  seismic  response  analysis  underestimates  the  block  forces 
caused  by  an  earthquake  the  magnitude  of  the  El  Centro 
earthquake.  The  one  degree  of  freedom  model  adequately 
predicts  the  vessel  rotational  response  to  seismic  loading  but 
is  ineffective  in  the  analysis  of  possible  failures  to  the 
vessel-drydock  system.  The  two  degree  of  freedom  model  does  a 
good  job  of  predicting  vessel  rotational  and  horizontal 
translation  response  to  seismic  loading  provided  the  response 
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due  to  vertical  acceleration  is  small,  which  is  the  case  for 
the  eleven  systems  analyzed.  Finally,  the  three  degree  of 
freedom  model  predicts  the  vessel  rotational,  horizontal,  and 
vertical  translation  responses  and  can  be  used  to  completely 
investigate  various  types  of  failure  to  the  vessel-drydock 
system  in  all  eleven  vessel-drydock  systems  examined.  The 
vessel  would  experience  side  pier  liftoff  failure  during  an 
earthquake  with  the  magnitude  of  the  El  Centro  earthquake  (ie. 
0.33  g's).  The  various  vessel-drydock  systems  would  not 
remain  intact  during  this  magnitude  of  earthquake.  Also  noted 
that  the  fourth  order  Runge-Kutta  numerical  scheme  should  be 
used  to  evaluate  system  response  because  of  its  simplicity, 
results  accurancy,  and  ease  of  introducing  terms  into  the 
Equation  of  motion. 

8 . 3  Recommendations  for  Further  Study 

The  seismic  response  of  drydocked  vessels  needs  further 
investigation  in  three  areas.  First,  a  three  degree  of 
freedom  model  which  will  allow  vessel  liftoff  should  be 
studied  so  that  other  failure  mode  will  occur.  Second,  a 
study  of  the  vessel-drydock  system  parameters  needs  to  be 
carried  out  to  explore  ways  in  which  drydocked  vessel  seismic 
response  can  be  decreased.  Third,  an  in  depth  study  of  the 
drydock  block  itself  in  order  to  provide  failure  modes  along 
with  stiffness  and  damping  characteristics  for  implementation 


into  mathematical  models. 
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APPENDIX  1 


LINEARIZED  THREE  DEGREE  OF  FREEDOM  SYSTEM  RESPONSE  USING 
MODAL  ANALYSIS  METHOD 

Integer  la, lb, ic, 1 ,mm,n, i job, iz, ier ,hul 1 , nsys , f laglO 
Integer  flagl, flag2, flag3, flag4. flagS, flag6, flag7, flag8 
real  a(6,6)  ,b (6,6)  ,beta,wk  (72)  , weight.h,  Ik,  gravity 
real  m(3,3) ,cx( 3,3) ,k(3,3) , sldearea .keelarea .plside 
real  bases ide, basekee 1 ,htside,htkeel , cr itl ,crlt2 ,cr lt3 
real  an(6,6)  ,ac(6000)  ,dtau,maxx,maxt,maxy,timex, timet 
real  r fl , r f2, r f 3,hfl ,hf2, hf 3, ampacc.mass, ampacmax 
real  kvs ,kvk,khs ,khk .base ,ht, counter 
real  timel , time 2 , time 3, time4, time 5 , time6 , time7, times 
real  x(6000) ,t(6000) ,y(60Q0) ,ri(7) ,si(7) ,pl(7) ,XSCL(6) 
real  bbb,ccc,wl2,wl,w22,w2,w32,w3,model ,mode3,capwldth 
real  mmxl .mmangl ,mmx3, mmang3,br , amp.plkeel , crlt4,ul,u2 
real  tlmey , mmmmml , mm.Timir.2  . mmmam3 , mmmmm4 
complex  al  fa  (6)  ,  z  (6 , 6)  ,  ad  (6 ,6)  ,bd  (6 , 6)  ,  a  a  (6 ,6)  ,bb  (6,6) 
complex  g(6)  ,v(6)  ,yy(6, 6000)  ,  ABC  (6)  ,  zt  (6,6)  ,betal 

CHARACTER *40  XLABEL , YLABEL , YYLAfiEL , YYYLABEL , DEC 


READ  IN  VESSEL  AND  DRYDOCK  DATA;  VESSEL  WEIGHT, KG, I (ABOUT  KEEL), 
TIME  INCREMENT  OF  DATA  POINTS, VERTICAL  STIFFNESS  OF  SIDE  AND 
KEEL  PIERS, HORIZONTAL  STIFFNESS  OF  SIDE  AND  KEEL  PIERS, 
CAVITATIONAL  CONSTANT, SIDE  BLOCK  BASE  AND  HEIGHT. 

KEEL  BLOCK  BASE  AND  HEIGHT, 

BLOCK-BLOCK  AND  BLOCK -HULL  FRICTION  COEFFICIENTS, 

SIDE  AND  KEEL  BLOCK'S  PROPORTIONAL  LIMIT, 

SIDE  PIER-VESSEL  CONTACT  AREA, KEEL  PIER-VESSEL  CONTACT  AREA, 

CAP  BLOCK  INCLINATION  ANGLE. 


read (44, *)  weight, h, Ik. dtau, kvs, kvk,khs,khk, gravity 
read (44,*)  baseside,  basekeel ,htslde,htkeel ,ul ,u2 
read (44, *)  br , amp.plslde.plkeel , sldearea , keel area, zet a 
read(44,*)  hull, nsys 

write  (6,*)  'do  you  want  response  plots?  (y  or  n) ' 

read (5, 15)  dec 

format  (a) 

do  10,1=1.6 

do  11, J=l,6 

a (1 , j) =0 .0 

b (1 , J) =0 .0 

continue 

continue 

do  12,1=1,3 

do  13, 3=1,3 

m (1, J) =0.0 

k (1 , j) =0 .0 

cx  (1 ,  j)  =0 .0 

continue 

continue 


CALCULATE  SYSTEM  PARAMETERS 


mass=veight/gravlty 

beta=asln (sqrt (br *  *2/ (4*h**2))) 

m  (1 , 1)  =mass 
m  (1 , 3)  =h*mass 
m(2, 2)  =mass 
m (3, 1) =mass*h 
m (3, 3) =Ik 

k (1 , 1)  =  (2*khs+khk) 
k (2 , 2) = (2*kvs+kvk) 

k  ( 3, 3)  =  ( (0 .5*kvs*  (br*  ‘2 .0) )  -  (welght'h) ) 

DETERMINE  NATURAL  FREQUENCIES  OF  SYSTEM 

bbb=-  (m  (1 , 1)  *k  (3,3)  +m(3,3)  *k  (l,l)-m(l,3)  *k  (3,l)-m(3,l)  *k  (1.3) ) 

/  (m  (1 , 1)  *m  (3, 3)  -m  (1 , 3)  *m(3,l)  ) 

ccc=  (k  (1, 1)  *k(3,3)-k(l,3)  *k(3,l))/(m(l,l)  *m  (3.  3) -m  (1 , 3)  *m(3,l)) 


NATURAL  FREQ.  MODE  #1 

wl2= (-bbb-sqrt (bbb* ‘2 .0-4*ccc) ) /2 
wl=sqrt  (wl2) 

NATURAL  FREQ.  MODE  #2 

v22=k(2,2)/m(2,2) 
w2=sqrt  (w22) 

NATURAL  FREQ.  MODE  #3 

w32= (-bbb+sqrt (bbb*  *  2 . 0-4*ccc) ) /2 
w3=sqrt  (w32) 

MODE  SHAPE  #1  &  #3 

model=  (m  (1 , 3)  *wl2-k  (1 , 3) )  /  (~m(l ,  1)  *wl2+k(l,l)) 

mode 3=  (m (1 , 3)  *w32-k  (1 , 3) )  /  (-m (1 , 1)  *w32+k(l,l) ) 

DETERMINE  Cll .C13.C31 ,C33 

mmxl=m(l,l)  +m  (1 . 3) /model 

mmangl=model*m (3, 1) +m (3, 3) 

mmx3=m  (1 , 1)  +m  (1 , 3)  /mode 3 

mmang3=mode3*m  (3, 1)  +m(3,  3) 

mnmmm  1 = 2  *  z  e  t  a  *  mmx  1  *  v  1 

mmmmm 2 = 2  * z e t a * mmx 3 * w 3 

mmmmm3=2*zeta*mmangl * vl 

mmmmm4=2*zeta*mmang3*v3 


cx  (1 , 3)  = (mmmmml-mmmmm2) / (l/model-l/mode3) 
cx  (1 , 1)  =mmmmml-  (cx  (1,3)  /model) 
cx(2,2)=2*zeta*m(2,2) *w2 


vvv.v 


continue 


READ  IN  ACCELERATION  DATA 
do  300 ,n=l , 5001 , 5 

read  (45,*)  ac  (n)  ,  ac  (n+1)  ,  ac  (n+2)  ,  ac  (n+3)  ,  ac  (n+4) 
continue 

ESTABLISH  FAILURE  CRITERIA  AND  FLAGS 

crltl=  min(ul,  (u2*cos  (beta) +sln  (beta) )  / 

(cos  (beta)  -u2*sin  (beta) ) ) 
crit2=mln  (ul,u2) 

cr It 3=  (0 .66*baseside-12 .0) /htside 

cr lt4=basekee 1/ (6  *htkeel ) 

ampacc=l .0 

counter=0 .0 

ampacmax=0 . 0 

continue 

flagl=0 

f lag2=0 

flag3=0 

f lag4=0 

f lag5=0 

flag6=0 

f lag7=0 

f lag8=0 

f lagl0=0 

do  50000  1=1,6 

ABC  (1)  =0.0 

continue 

maxx=0 . 0 

maxt=0 .0 

maxy=0 . 0 

mm=0 

SOLVE  FOR  Y.THE  COLUMN  MATRIX  WHICH  IS  THE  SOLUTION  IN 

THE  IMAGINARY  COORDINATE  SYSTEM  WHERE  THE  EQUATIONS  OF  MOTION 

ARE  DECOUPLED 

do  301,1=3,3501,2 
do  302,1=1,6 

g  (1)  =  (ABC  (1)  +  ac  (1-2)  /2 . 54)  *exp(v(i)  *2*dtau) 

+4.0*  (ac  (1-1)  /2 .54)  *exp(v(i)  *dtau)  +ac  (l)/2.54 
ABC  (1)  =g  (1) 

yy  (1, 1)  =  ( (zt  (1,4)  +amp*zt  (1 , 5)  +zt  (1 ,6)  *h)  *nass/ad  (1 , 1J )  *g  (1) 

*ampacc*-dtau/3 

continue 

mm=mm+l 

USING  THE  MODAL  MATRIX  (  4th, 5th,  AND  6th  ROWS  ),  OBTAIN 
VALUES  FOR  TRANSLATIONS  AND  ROTATION 

do  303,  1=1,6 

si(l)=z(4,l)*yy(i,l)+sl(l-l) 
continue 
x (nun)  =sl  (6) 

If  (abs  (x  (miri) )  . gt . abs  (ir.axx) )  then 
tlmex=dtau* (1-1) 
n:axx=x  (irun) 

endl  f 

do  304.  1=1,6 


n  n 


ri(i)=z  (6,l)*yy(i,l)+rl(l-L) 
continue 
t  (mm)  =ri  (6) 

If  (abs  (t  (mm) )  .gt .  abs  (maxt) )  then 
timet=dtau* (1-1) 
maxt=t (mm) 

endif 

do  305,  1=1,6 

pl(l)=z(5,l) *yy(l,l)+pi(l-l) 
continue 
y  (mm)  =pl  (6) 

If  (abs  (y  (mm) )  .  gt .  abs  (maxy) )  then 
timey=dtau* (1-1) 
maxy=y  (mm) 

endif 

CALCULATE  VERTICAL  AND  HORIZONTAL  FORCES  CAUSED  BY  VESSEL, 

TEST  FOR  FAILURE 

CALCULATE  FORCES  ON  SIDE/KEEL  BLOCKS 
r fl=kvs*  ( (veight/k  (2,2))  -y  (mm)  -  (br/2)  *t  (mm)  ) 
r  f  2=kvs*  (  (veight/k  (2 , 2)  )  -y  (mm)  +  (br/2)  *t  (mm)  ) 
rf3=kvk* ( (veight/k (2, 2) ) -y (mm) ) 
hfl=khs*  (x  (mm) ) 
hf2=khs*  (x  (mm) ) 
hf 3=khk*x (mm) 

TEST  FOR  SIDE  BLOCK  SLIDING 

if  (flagl.eq.l)  then 
go  to  400 

else  if  (hfl .  It  .0 .0  .and.r  f  1  .gt  .0 .0  .and. abs  (hfl/r  fl)  .gt  .critl)  then 
timel=  dtau*  (1-1) 
f lagl=l 

else  if  (hf  2  .  gt  .0 .0  .and.r  f2  .gt  .0 .0  .and. abs  (hf2/rf2)  .gt  .critl)  then 
timel=dtau*  (1-1) 
f lagl=l 

endif 
xl=x  (mm) 
yl=y (mm) 
tl=t  (mm) 
continue 

TEST  FOR  KEEL  BLOCK  SLIDING 

if  (flag2.eq.l)  then 
go  to  410 

else  if  (r  f3  .gt  .0 .0  .  and.  abs  (hf  3/r  f  3)  .gt  .crit2)  then 
time2=dtau* (1-1) 
f  lag2=l 

endif 
x2=x (mm) 
y2=y (mm) 
t2=t (mm) 
continue 

TEST  FOR  SIDE  BLOCK  OVERTURNING 

if  (flag3.eq.l)  then 
go  to  420 

else  if  (hf 1 . It .0 . 0 . and .r fl . gt .0 .0 . and. abs (hfl/r fl) . gt .crlt3)  then 
tlme3=  dtau* (1-1) 
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f lag3=l 

else  If  (hf2  .gt  .0 .0  .and  .r  f2  .gt  .0 .0  .  and.abs  (hf2/rf2)  .gt  .crlt3)  then 
tlme3=dtau* (1-1) 
f lag3=l 

endif 
x3=x (mm) 
y3=y  (mm) 
t3=t (mm) 
continue 

TEST  FOR  KEEL  BLOCK  OVERTURNING 

If  (flag4.eq.l)  then 
go  to  430 

else  If  (r  f  3 .  gt  .0 .0  .and  .abs  (hf  3/r  f 3)  .gt  .crlt4)  then 
tlme4=dtau* (1-1) 
f lag4=l 

endif 
x4=x (mm) 
y 4=y (mm ) 
t4=t (mm) 
continue 

TEST  FOR  SIDE  BLOCK  LIFTOFF 

If  (flagS. eq.l)  then 
go  to  440 

else  If  (rfl.lt. 0.0  .or.  rf2.lt. 0.0)  then 
time5=dtau* (1-1) 
f lag5=l 

endif 
x5=x (mm) 
yS=y  (mm) 
tS=t (mm) 
continue 

TEST  FOR  KEEL  BLOCK  LIFTOFF 

If  (flag6.eq.l)  then 
go  to  450 

else  If  (rf3.lt. 0.0)  then 
tlme6=dtau* (1-1) 
f lag6=l 

end!  f 
x6=x (mm) 
y6=y  (mm.) 
t6=t (mm) 
continue 

TEST  FOR  SIDE  BLOCK  CRUSHING 

If  (flagV.eq.l)  then 
go  to  460 

else  If  (rfl.gt.0.0  .and.  (r f 1/sldearea) . gt .plside)  then 
f lag7=l 

time7=dtau* (1-1) 

else  if  (rfZ.gt.0.0  .and.  (rf 2/sldearea) . gt .plside)  then 
flag7=l 

time7=dtau* (1-1) 


endi  f 
x7=x  (mis) 
y7=y  (mm) 
t7=t (mm) 
continue 

TEST  FOR  KEEL  BLOCK  CRUSHING 

If  (flag8.eq.l)  then 
go  to  470 

else  if  (rf3.gt.0.0  .and.  (r f 3/keelarea) .gt .plkeel)  then 
flag8=l 

time8=dtau*  (1-1) 

end  If 
x8=x  (mm) 
y8=y (mm) 
t8=t (mm) 
continue 

continue 


PLOT  RESULTS 
go  to  999 
continue 

if  (DEC.EQ. 'N' )  THEN 

write (6.*)  'I  am  finishing.' 

GO  TO  998 

endi  f 

write (6, *)  'I  am  plotting.’ 

XSCL (1) =0 . 0 
XSCL  (2)  =30 . 0 
XLABEL=’TIME  IN  SECONDS' 

YLABEL= 'ROTATION  IN  RADIANS' 

YYLABEL= 'RELATIVE  HORIZONTIAL  DISPLACEMENT  IN  INCHES' 
YYYLABEL= 'RELATIVE  VERTICAL  DISPLACEMENT  IN  INCHES' 

CALL  QPICTR  (X . 1 . 1500 .QXSCL (XSCL) ,QISCL (-2) -QXLAB (XLABEL) 
QYLAB (YYLABEL)  ,QLABEL(4)  ) 

CALL  QPICTR  (T. 1 , 1500 .QXSCL (XSCL) ,QISCL (-2) ,QXLAB (XLABEL) 
QYLAB (YLABEL) .QLABEL (4)  ) 

CALL  QFICTR  (Y. 1 , 1500 .QXSCL (XSCL) ,QISCL (-2) .QXLAB (XLABEL) 
QYLAB (YYYLABEL) .QLABEL (4) ) 
go  to  200C0 
CONTINUE 


i  f  (ampacc .  eq .  1 . 0)  then 
write  (46, 4000)  r.sys 

format (lx,/, 28x, '**•*  System  ’ , 1 2 , lx, ' *  *  *  * ' ) 
write  (46 , 4050)  hull 

format (lx, /, 30x, ' *  *  Hull  ' , 1 3. lx, ’ *  * ' ) 
write (46 , 4100) 

format (lx. //. 28x , ' *  Ship  Parameters  *') 
write (46 . 4150)  , 

format(lx,/,5x, 'Weight' ,8x, 'Momt.it  of  Inertia', 9x,  K.G.  ) 
write  (46,4200)  weight. Ik, h 

format  (lx . f 9 . 1 , lx . ' kips ' . lx. f 11 . 1 , lx, 'kips-ln-sec2 ' . 

+3x, f6 . 1 , IX, 'ins') 
write (46, 4250) 

format (lx.//, 26x. ’ *  Drydock  Parameter  *’) 
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write  (46 , 4300) 

format (lx, /, lx, ' Side  Block  Height *, 3x,  1  Side  Block  Width’, 

►  3x, ’Keel  Block  Height 3x, 'Keel  Block  Width’) 

write  (46 , 4350)  htside.baseslde.htkeel , basekeel 

format (2x, f6 . 1 , lx. 1  ins ’ , llx, f6 . 1 , lx, *  ins ' , llx, f6 . 1 , lx, ’ ins ’ , 

►  9x, f6 . 1, lx, 1  Ins ’ ) 

write  (46,4400) 

format (lx, /, lx, ’ Side-to-Side  Pier  Distance’) 

write (46 , 4450)  br 

format (lx, t7 , f6 . 1 . lx, ’ins’) 

write  (46,4470) 

format (lx,/, *  Total  Side  Pier  Contact  Area’ 

►,3x, 'Total  Keel  Pier  Contact  Area’) 
write  (46 , 4475)  sidearea.keelarea 
format (Ix,8x,fll.l,lx,'ln2’,14x,fll.l,lx,'ln2’) 
write  (46,4500) 

format (lx, /, lx, ' Block-on -Block  Friction  Coef f ' , 3x, 'Hul i-on-Biock 

►  Friction  Coeff') 

write  (46 , 4550)  ul,u2 
format (lx, lOx, f7 .3, 23x, f7. 3) 
write  (46,4600) 

format (lx,/, lx, 'Side  Pier  Proportional  Limit' , 3x, 'Keel  Pier' 
Proportional  Limit') 
write (46 , 4650)  plside,plside 

format (lx, lOx, f7 . 3, lx, 'kips/in2 ' 15x, f7 . 3, lx, 'kips/ in 2 ' ) 
write (46 , 4700) 

format (lx,/, lx, ' Side  Pier  Vertical  Sti f fness ' , 3x, ' Side  Pier', 

►’  Horizontal  Stiffness’) 
write(46,4750)  kvs.khs 

format (lx, 3x, f 11 . 1 , lx, 'kips/in' , llx, f 11 . 1 , lx, 'klps/ln ' ) 
write (46, 4775) 

format (lx,/, lx, 'Keel  Pier  Vertical  Stif fness ’, 3x, 

►'Keel  Pier  Horizontal  Stiffness') 
write (46, 4780)  kvk.khk 

format (lx, 3x, f 11 . 1 , lx, 'kips/in' , llx, f 11 . 1, lx, 'klps/in ' ) 
write  (46,4800) 

format (lx,//, 20x, ' *  System  Parameters  and  Inputs  *') 
write  (46,4850) 

format (lx,/, lx, ' Input  Forcing  Function  is  Horizontal  Component', 
►'  of  the  1946  El  Centro') 
write  (46 , 4875) 

format (lx, 20x, '  Earthquake  Acceleration  Time  History.') 
write  (46,4995) 

format (lx,/, lx, 'Vertical/Horizontal  Ground  Acceleration  Ratio’ 
►,3x,'Data  Time  Increment') 
write  (46 , 4990)  amp.dtau 
format (lx,10x, f6. 3,tS5, f6.3,lX, 'sec') 
write (46, 4900) 

format (lx,/, lx, 'Gravitational  Constant' , 3x, 'Percent  System 

►  Damping') 

write (46 , 4950)  gravity, zeta*100. 

format (lx,7x, f6.2,ix, ' in/sec 2 ' , lOx, f6.2,lx, '%') 

write (46, 5000) 

format (lx , /, 25x , ' Mass  Matrix ', /) 
do  5100  1=1,3 

write  (46 , 5050)  m  (i  ,  1)  ,m  (i  ,  2)  ,m  (i  ,  3) 
format (lx, flS.4. Sx. flS.4,5x, fl5.4) 
continue 
write  (46 , 5200) 
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5200 


5250 

5300 


5400 


5450 

5S00 


5700 

6000 


6001 


6002 


6500 


format (lx. /, 25x, 'Damping  Matrix1 ,/) 
do  5300  1=1,3 

write (46, 5250)  cx(l,l)  ,cx(l,2)  ,cx(l,3) 
format (lx, fl5.4,5x, fl5.4, 5x, f 15 .4) 
continue 
write (46,5400) 

format (lx./, 25x, ’Stiffness  Matrix ’ ,/) 
do  5500  1=1.3 

write (46,5450)  k(i.l)  ,k(l,2)  ,k(l,3) 

format (lx, f 15 . 4, 5x, f 15 . 4, 5x, fl5.4) 

continue 

write (46,5700) 

format (lx,//) 

WRITE (46,6000) 

FORMAT (IX, 'Undamped  Natural  Frequencies ', t35, ’Mode  #l’,t50, 

+'Mode  #2 ’ , t65, ‘Mode  #3') 
write  (46,6001)  wl,w2,w3 

format (lx, t31 , f7 . 3, lx, ' rad/sec ' , t46 , f7 . 3, lx, ’rad/sec ’ , t62 , f7 . 3, 

+ ’  rad/sec ’ ) 

WRITE (46,6002) 

FORMAT (IX, 'Damped  Natural  Frequencies ', t35, 'Mode  #l',t50, 

♦'Mode  #2' ,t65, 'Mode  #3') 

WRITE (46,6500)  wl‘sqrt (l-zeta**2) ,w2*sqrt (1-zeta*  *  2) , w3*sqrt (l-zeta**2) 
format (lx, t31 , f7 . 3, lx, 'rad/sec ' ,t46, f7 . 3, lx, 'rad/sec' ,t62, f7 . 3, 

+  '  rad/sec') 
endif 


10500 


write (46, 10500)  ampacc*100 

format (lx,///, lx, 'For  Earthquake  Acceleration  of  ',f6.2,'  % 
'of  the  El  Centro',/) 


25000 


write (46, 25000) 

format (lx, 'Maxlmums/Fallures ', t26, 'X  (ins) ' , t36, ' Y  (ins)',t51, 

♦  'Theta  (rads)  '  ,t65, 'Time  (sec)') 
write (46, 25001) 

25001  format  (lx, ' - ’  ,t25,  ’ - '  ,t35,  ’ - ’  ,t50, 

+  ' - ’  ,t64,  ' - ') 

write  (46,310)  maxx.tlmex 
format  (lx, *  Maximum  X* , t25, f9 .6, t65, f5 . 2) 
write  (46,311)  maxy.tlmey 
format  (lx,1  Maximum  Y' , t35, f9 .6 , t65, f5 . 2) 
write  (46,312)  maxt, timet 

format  (lx,'  Maximum  Rotation ', t50 , f9 .6, t65, f5 . 2) 


310 

311 

312 


313 


If  (flagl.eq.l)  then 
f lagl0=f Iagl0+1 

write  (46,313)  xl ,yl , tl , tlmel 

format  (lx, 'Side  block  sliding'  , t25, f9 . 6, t35, f9 .6 , t50 , f 9 .6 , 
+t65, f5 . 2) 


endif 


314 


If  (flag2.eq.l)  then 
f lagl0=f laglO+l 

write  (46,314)  x2 ,y2 , t2 , time2 

format  (lx, 'Keel  block  sliding’  , t25, f9 .6, t3S, f9 .6, t50, f9 .6, 
+t65, f 5 . 2 ) 
endif 


if  (flag3.eq.l)  then 
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f laglO=f laglO+1 

write  (46,315)  x3,y3.t3,time3 

315  format  (lx, 'Side  block  overturning’  , t25, f9 .6, t35. f9 .6, t50, f9 .6 , 
+t65, f5. 2) 

endif 

if  (flag4.eq.l)  then 
flaglO=flaglO+l 

write  (46,316)  x4,y4,t4, time4 

316  format  (lx,'Keel  block  overturning'  , t25. f9 .6, t35, f9 .6 , t50 , f9 .6 , 
+t65, f5 .2) 

endif 

if  (flag5.eq.l)  then 
f laglO=f laglO+1 

write  (46,317)  x5,y5, t5, time5 

317  format  (lx, ' Side  block  liftoff'  , t25, f9 .6, t35, f9 .6, tSO , f9 .6 , 
+t65,  f5 . 2) 

endif 


* 


if  (flag6.eq.l)  then 
f  laglO=f  laglO-^1 

write  (46,318)  x6, y6, t6, time6 

318  format  (lx, 'Keel  block  liftoff '  . t25. f9.6. t35. f9 .6. tSO, f9 .6, 
+t65, f 5  -  2 ) 
endi  f 


M 


if  (flag7.eq.l)  then 
f laglO=f laglO+1 

write  (46,319)  x7,y7,t7,time7 

319  format  (lx, 'Side  block  crushing’  , t2S, f9 .6, t35, f9 .6, t50 , f9 .6 , 

+t65, fS. 2) 
endif 


£ 


.** 


if  (flag8.eq.l)  then 
f laglO=f laglO+1 

write  (46,320)  x8 ,y8, t8. time8 

320  format  (lx, ’ Keel  block  crushing’  , t25, f9 .6, t35, f9 .6, t50, f9 .6. 

+t65,  f5 . 2) 
endif 

if (flaglO.le.O)  then 
write  (46, 11000) 

11000  format (lx,/, lx, ’No  failures  occurred.’) 

1 f (counter ,eq . 1 . 0  .and.  f laglO . la .0)  then 

go  to  60000 

endif 

if (counter .eq .0 .0)  then 
ampacmax=arapacc 
ampacc=ampacc+ . 1 
counter=l .0 

wrlte(6,*)  ’I  am  in  the  secondary  looping  stage.’ 

endl  f 

endif 

i f (ampacc . le . ampacmax)  go  to  20000 

if (counter .eq. 1 .0)  then 

ampacc=ampacc- .01 

else  1 f (counter .eq . 0 .0)  then 

ampacc=ampacc- . 1 

endi  f 


M 


Vv< 
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go  to  10000 

continue 

stop 

end 

subroutine  mult  (a,b,cyy) 
real  b(6,6) 

complex  a(6,6)  ,cyy(6,6)  ,d(6) 

do  100, 1*1,6 

do  200, J=l,6 

do  300 ,k=l , 6 

d  (k)  =a  (1  ,k)  *b(k,J) 

continue 

cyy  (1,  J)  =d  (1)  +d  (2)  -d  (3)  -d  (4)  *d  (5)  +d  (6) 

continue 

continue 

return 

end 

subroutine  multc  (a.b.cyx) 

complex  a(6,6) ,b(6,6) ,cyx(6,6) ,d(6) 

do  100,  1=1,6 

do  200,  j=l , 6 

do  300,  k=l ,6 

d(k)=a(l,k)  *b  (k,  j) 

continue 

cyx  (1,  J)  =d  (1)  -*-d  ( 2 )  +d(3)->-d(4)  +d(5)+d(6) 
continue 
continue 
return 
:nd 


APPENDIX  2 


NON-LINEAR  THREE  DEGREE  OF  FREEDOM  SYSTEM  RESPONSE 
USING  FOURTH  ORDER  RUNGE-KUTTA  METHOD 

Integer  NN, 1 ,mm,n,hul 1 .nsys , flaglO 

Integer  flagl, flag2,  flag3, flag4, flag5, flag6, flag7,  flag8 
real  beta , weight ,h, Ik, gravity 

real  m(3,3)  , cx  (3, 3)  ,k  (3, 3)  ,  si  dearea ,  keel  area  , plside 

real  baseside,basekeel,htside,htkeel,critl,crlt2,crit3 

real  ac(6000)  ,dtau,maxx,maxt,maxy,timex, timet 

real  rfl,rf2,rf3,hfl ,hf 2 ,hf 3, ampacc .mass , ampacmax 

rea 1  kvs , kvk , khs . khk , base , ht , counter 

real  timel , t ime2, time 3, time4, time5, t ime6 , time7, time8 

real  x(6000)  ,t(6000)  ,y(6000)  ,XSCL(6) 

real  bbb, ccc, wl2 ,wl ,w22, w2 , w32 . w 3, model ,mode3,capwidth 
real  mmxl /ininangl ,nimx3 ,mmangj,br , amp.plkeel , crit4,ul,u2 
real  timey,mninmml,nramnm2,i7mimmni3,mmnmm4 
real  R,S,TAU.A(5)  ,B(5)  ,0(5)  ,D(5)  ,E(5)  ,F(5)  ,C(5)  ,HH(5) 
CHARACTER *40  XLABEL, YLABEL, YYLABEL, YYYLABEL,DEC 


READ  IN  VESSEL  AND  DRYDOCK  DATA:  VESSEL  WEIGHT, KG, I (ABOUT  KEEL) 
TIME  INCREMENT  OF  DATA  POINTS, VERTICAL  STIFFNESS  OF  SIDE  AND 
KEEL  PIERS, HORIZONTAL  STIFFNESS  OF  SIDE  AND  KEEL  PIERS, 

GAV I TAT I ONAL  CONSTANT. SIDE  BLOCK  BASE  AND  HEIGHT, 

KEEL  BLOCK  BASE  AND  HEIGHT , 

BLOCK-BLOCK  AND  BLOCK-HULL  FRICTION  COEFFICIENTS, 

SIDE  AND  KEEL  BLOCK'S  PROPORTIONAL  LIMIT, 

SIDE  PIER-VESSEL  CONTACT  AREA, KEEL  PIER-VESSEL  CONTACT  AREA, 

CAP  BLOCK  INCLINATION  ANGLE. 

read  (44, *)  welght.h, Ik.dtau, kvs, kvk, khs, khk, gravity 
read(44,*)  baseslde,  basekeel , htslde.htkeel ,ul ,u2 
read (44, *)  br , arap.plslde ,plkeel , sidearea, keel area, zeta 
read (44,*)  hull, nsys 

write  (6,*)  'do  you  want  response  plots?  (y  or  n) ' 

read  (5, 15)  dec 

format  (a) 

do  12,1=1,3 

do  13, 3=1 , 3 

m  ( 1 ,  J )  =0 . 0 

k ( i , J ) =0 . 0 

cx(l, J)=0.0 

continue 

continue 


CALCULATE  SYSTEM  PARAMETERS 


mass  =we 1 ght/gr  a v 1 ty 
beta=asln (sqrt (br *  *2/ (4*h*  *  2) ) ) 


.  ,\ 


n  o 


read  (45,*)  ac  (n)  ,  ac  (n+1)  ,  ac  (n->2)  ,  ac  (n*  3)  ,  ac  (n+4) 
300  continue 

C  ESTABLISH  FAILURE  CRITERIA  AND  FLAGS 

critl=  min  (ul ,  (u2*cos  (beta)  ♦sin  (beta) )  / 

♦  (cos (beta) -u2*sin (beta) ) ) 
crit2=min  (ul,u2) 

crit3=  (0 .66*baseside-12 .0) /htside 
crit4=basekeel/(6*htkeel) 
ampacc=l .0 
counter =0 .0 
ampacmax=0 .0 
10000  continue 
flagl=0 
f lag2=0 
flag3=0 
f lag4=0 
f lag5=0 
f lag6=0 
flag7=0 
f lag8=0 
f lagl0=0 
maxx=0 . 0 
maxt=0 .0 
maxy=0 .0 
mm=0 

x(l)=0.0 

y(i)=o.o 

t(l)=0.0 
R=0 .0 
S=0.0 
TAU=0 . 0 


IMPLEMENTATION  OF  EQUATIONS  OF  MOTION  INTO  THE 
RUNGE-KUTTA  FOHMULUS 

do  301,1=1,6000 
DO  3000,11=0,4 
A (1 1) =0 .0 
B ( 1 1 ) =0 . 0 
C(11)=0.0 
D (11) =0 .0 
E ( 1 1 ) =0 . 0 
F (1 1) =0 .0 
G ( 1 1 ) =0 . 0 
HH(ll)  =0 .0 
3000  CONTINUE 
mm=mm*  1 

DO  302,  NN=1 , 4 
IF (NN.EQ. 1)  THEN 
FF=0 . 0 

ELSE  IF  (NN.EQ. 2  .OR.  NN.EQ. 3)  THEN 
FF=0 . 5 

ELSE  IF  (NN.EQ.  4)  THEN 
FF=1 .0 
ENDIF 

A(NN)  =dtau*  (R+FF*D(NN-1) ) 

B (NN) =dtau* (S*FF‘E (NN-1) ) 


t.'-r.r 
/  V  > 


.v.sV 

■  V  v 
v'V, 

•  v  »* 1 
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C (NN) =dtau* (TAU+FF*F (NN-1) ) 

D(NN)=dtau*  ( (-cx  (2 , 2) /m (2 , 2) )  *  (R+FF*D  (NN-1) )  -  (k  (2 , 2)  /m  (2 , 2) ) 

+  *  (y  (mm)  +FF*A  (NN-1) )  -amp*ampacc*ac  (1)  /2 .54) 

G  (NN)  =dtau»  ( (-cx  (1 . 1)  /m  (1 . 1) )  *  (S+FF*E  (NN-1) )  -  (cx  (1 , 3)  /m  (1 , 1) ) 

+  *  (TAU+FF*F  (NN-1) )  -  (k  (1 , 1)  /m  (1 , 1) )  *  (x  (mm) +FF*B (NN-1) ) -ampacc*ac (1) / 
♦2.54) 

HH (NN)  =dtau*  ( (-cx  (3, 3)/m(3. 3) )  *  (TAU+FF*F  (NN-1) )  -  (cx  (3, 1)  /m  (3. 3) ) 
+  *  (S+FF*E  (NN-1) )  -  (k  (3, 3)  /m  (3 , 3) )  *  (t(mm)  +FF‘C  (NN-1) )  ♦  (m  (3, 1) /m  (3, 3) ) 
+  *  (  (-cx  (2 . 2)  /m  (2 , 2) )  *  (R+FF*D  (NN-1)  )  -  (k  (2 , 2)  /m  (2 , 2)  )  *  (y  (mm)  +FF*A  (NN- 
♦  1) )  )  *  (t  (mm)  +FF*C  (NN-1) )  -  (m  (3, 1)  /m  (3. 3)  )  *ampacc*ac  (l)/2.54) 


E  (NN)  =  (m  (1 , 1)  *m  (3,  3)  *C  (NN)  -m  (1 , 3)  *m(3,3)  *HH  (NN)  )  / 
+  (m  (3, 3)  *m  (1 . 1)  -m  (1 , 3)  *m(3. 1) ) 

F  (NN)  =  (HH(NN)  -  (m  (3. 1)  /m  (3,  3) )  *E  (NN) ) 

302  continue 


DETERMINING  SYSTEM  RESPONSE 
y  (mm+1)  =y  (mm)  +  (A(l)  +2*A(2)  ♦2*A(3)  ♦A(4))/6 
x  (mm+1)  =x  (mm)  +  (B(l)+2*B(2)+2*B(3)+B(4))/6 
t  (mm+1)  =t  (mm)  +  (C(l)+2*C(2)+2*C(3)+C(4))/6 
R=R+  (D (1)  +2*D (2)  +2*D(3)  +D(4)  )/6 
S=S+(E(l)+2*E(2)+2*E(3)+E(4))/6 
TAU=TAU+  (F  (1)  +2*F  (2)  +2*F  (3)  +F  (4)  )  /6 


MAXIMUM  VALUES  FOR  TRANSLATIONS  AND  ROTATION 

If  (abs  (x  (mm) )  .  gt .  abs  (maxx) )  then 
tlmex=dtau* (1-1) 
maxx=x (mm) 

endl  f 

If  (abs  (t  (mm)  )  .gt .  abs  (maxt)  )  then 
tlmet=dtau* (1-1) 
maxt=t (mm) 

endl  f 

If  (abs  (y  (mm) )  . gt .  abs  (maxy) )  then 
tlmey=dtau* (1-1) 
maxy=y (mm) 

endl  f 

CALCULATE  VERTICAL  AND  HORIZONTAL  FORCES  CAUSED  BY  VESSEL, 
TEST  FOR  FAILURE 

CALCULATE  FORCES  ON  SIDE/KEEL  BLOCKS 
r  f  l=kvs  *((weight/k(2,2))-y(mm)-  (hr/2)  *t  (mm) ) 
r  f2=kvs*  ( (veight/k  (2, 2) )  -y  (mm)  +  (br/2)  *t  (mm) ) 
r f3=kvk* ( (welght/k (2,2)) -y (mm) ) 
hf  l=khs*  (x  (mm)  ) 
hf  2=khs*  (x  (mm) ) 
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hf3=khk*x  (mm) 

TEST  FOR  SIDE  BLOCK  SLIDING 


» 

I  • 

I  400 

I  C 


410 

C 


i 


420 

C 


I 


430 


C 


If  (flagl.eq.l)  then 
go  to  400 

else  If  (hfl . It .0 .0 .and .rfl .gt .0 .0 .and. abs (hfl/r fl) .gt .crltl)  then 
tlmel=  dtau* (1-1) 
f lagl=l 

else  if  (hf2 .gt .0 .0 .and.r f2 .gt .0 .0 .and. abs (hf2/r f 2) .gt .crltl)  then 
tlmel=dtau* (1-1) 
f lagl=l 

endl  f 
xl=x  (mm) 
yl=y (mm) 
tl=t  (mm) 
continue 

TEST  FOR  KEEL  BLOCK  SLIDING 

If  (flag2.eq.l)  then 
go  to  410 

else  If  (r  f3 .  gt . 0 . 0 .  and  .  abs  (hf 3/r  f 3)  .gt  .crit2)  then 
tlme2=dtau* (1-1) 
f lag2=l 

endl 

x2=x (mm) 
y2=y (mm) 
t2=t (mm) 
continue 

TEST  FOR  SIDE  BLOCK  OVERTURNING 

If  (flag3.eq.l)  then 
go  to  420 

else  If  (hfl .  It  .0 .0  .  and  .rfl  .gt.0 .0  .and. abs  (hfl/r  fl)  .gt  .crlt3)  then 
tlme3=  dtau* (1-1) 
f lag3=l 

else  If  (hf  2  .gt  .0 .0  .  and.r  f2  .gt  .0 .0  .and.abs  (hf2/r  f2)  .gt  .crlt3)  then 
tlme3=dtau* (1-1) 
f lag3=l 

endlf 
x3=x (mm) 
y3=y  (mm) 
t3=t (mm) 
continue 

TEST  FOR  KEEL  BLOCK  OVERTURNING 

If  (flag4.eq.l)  then 
go  to  430 

else  If  (r  f 3 . gt . 0 .0 .  and .  abs  (hf 3/r  f 3)  .  gt . crit4)  then 
time4=dtau*  (1-1) 
f lag4=l 

endl  f 
x4=x (mm) 
y4=y (mm) 
t4=t (mm) 
continue 

TEST  FOR  SIDE  BLOCK  LIFTOFF 
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if  (flag5.eq.l)  then 
go  to  440 

else  If  (rfl.lt. 0.0  .or.  rf2.lt. 0.0)  then 
time5=dtau* (1-1) 
flag5=l 

endlf 
x5=x (mm) 
y5=y  (mm) 
t5=t  (mm) 

440  continue 

C  TEST  FOR  KEEL  BLOCK  LIFTOFF 

if  (flag6.eq.l)  then 
go  to  4S0 

else  if  (rf3.lt. 0.0)  then 
tlme6=dtau*  (1-1) 
f lag6=l 

endif 
x6=x (mm) 
y6=y (mm) 
t6=t (mm) 

450  continue 

C  TEST  FOR  SIDE  BLOCK  CRUSHINC 

if  (flag7.eq.l)  then 
go  to  460 

else  if  (rfl.gt.0.0  .and.  (r f 1/sidearea) . gt .plside)  then 
f lag7=l 

time7=dtau*  (1-1) 

else  if  (rf2.gt.0.0  .and.  (rf2/sidearea) .gt. plside)  then 
f lag7=l 

tlme7=dtau*  (1-1) 

endif 
x7=x  (mm) 
y7=y  (mm) 
t7=t (mm) 

460  continue 

C  TEST  FOR  KELL  BLOCK  CRUSHING 

if  (flag8.eq.l)  then 
go  to  470 

else  if  (rf3.gt.0.0  .and.  (r f 3/keel area) . gt .plkeel)  then 
f lag8=l 

tlme8=dtau* (1-1) 

end!  f 
x8=x (mm) 
y8=y  (mm) 
t8=t (mm) 

470  continue 

301  continue 


C  PLOT  RESULTS 

go  to  999 
60000  continue 


if  (DEC.EQ. 'N')  THEN 

write  (6,*)  'X  am  finishing.' 

CO  TO  998 

endlf 

wrlte(6,*)  'I  am  plotting.' 

XSCL (1) =0 .0 

XSCL (2) =30 .0 

XLABEL= 'TIME  IN  SECONDS' 

YLABEL= 'ROTATION  IN  RADIANS' 

YYLABEL=' RELATIVE  HORIZONTIAL  DISPLACEMENT  IN  INCHES' 

YYYLABEL= 'RELATIVE  VERTICAL  DISPLACEMENT  IN  INCHES' 

CALL  QPICTR  (X, 1 , 3000 .QXSCL (XSCL) .QISCL (-2) , QXLAB (XLABEL) , 

+  QYLAB (YYLABEL) .QLABEL (4)  ) 

CALL  QPICTR  (T, 1 , 3000 , QXSCL (XSCL) ,QISCL(-2) , QXLAB (XLABEL) , 

+  QYLAB (YLABEL) .QLABEL (4)) 

CALL  QPICTR  (Y. 1 , 3000 .QXSCL (XSCL) ,QISCL (-2) . QXLAB (XLABEL) . 

*  QYLAB (YYYLABEL) .QLABEL (4)  ) 
go  to  20000 
CONTINUE 

1 f (ampacc . eq. 1 .0)  then 
write (46,4000)  nsys 

format (lx, /, 28x, '*** *  System  ', 12, lx, '****’ ) 
write  (46 , 4050)  hull 

format (lx,/, 30x, '**  Hull  ' , 13, lx, ' *  * ’ ) 
write  (46.4100) 

format (lx,//, 28x, ' *  Ship  Parameters  *') 
write (46, 4150) 

format (lx,/, 5x, 'Weight' ,8x, 'Moment  of  Inertia' ,9x, 'K.G. ' ) 
write (46,4200)  weight, Ik, h 

format (lx, f9 . 1, lx. 'kips ' , lx, f 11 . 1 , lx, ’kips-ln-sec2 ’ , 

+3x, f6 . 1 > IX, 'Ins'} 
write  (46, 4250) 

format (lx,//, 26x. ' *  Drydock  Parameter  *') 
write (46, 4300) 

format (lx, /, lx, ' Side  Block  Height ', 3x, ' Side  Block  Width’, 
+3x,'Keel  Block  Height' , 3x, 'Keel  Block  Width’) 
write (46 , 4350)  ht side .bases lde.htkeel , basekeel 
format (2x, f6 . 1 , lx, ' Ins ' , llx, f6 . 1 , lx, ' Ins ' , llx, f6 . 1 , lx, ' Ins ' , 

+9x, f6 . 1 , lx, ' ins ’ ) 
write (46,4400) 

format (lx,/, lx, ' Side-to-Slde  Pier  Distance') 

write (46.4450)  br 

format (lx, t7 , f6 . 1 , lx, ' Ins ' ) 

write (46. 4470) 

format (lx,/, '  Total  Side  Pier  Contact  Area' 

+,3x, 'Total  Keel  Pier  Contact  Area') 
write  (46 , 4475)  sidearea.keelarea 
format (lx, 8x , f 11 . 1 , lx, ' ln2 ' , 14x, fll . 1 , lx, ’ ln2 ' ) 
write (46,4500) 

format (lx,/. lx, 'Block-on-Block  Friction  Coeff ' , 3x, 'Hul 1-on-Block 
+  Friction  Coeff) 

write  (46 , 4550)  ul,u2 
format ( lx , lOx , f 7 . 3 , 2  3x . f 7 . 3) 
write (46, 4600) 

format (lx,/. lx, ' Side  Pier  Proportional  Limit ', 3x, ' Keel  Pier' 

+,'  Proportional  Limit') 

write (46, 4650)  plslde.plslde 

format (lx,10x,f7.3,lx, ' klps/in2 ' 15x, f7 . 3, lx, 'klps/ln2 ’ ) 


write  (46,4700) 

format (lx,/, lx, ' Side  Pier  Vertical  Stif fness ' , 3x, ' Side  Pier’, 

+  '  Horizontal  Stiffness') 
write  (46 , 4750)  kvs,khs 

format (lx, 3x, fll . 1 , lx, 'kips/in ' , llx, fll . 1 , lx, 'kips/in ’ ) 
write  (46,4775) 

format (lx,/, lx, 'Keel  Pier  Vertical  Stif fness ', 3x, 

+'Keel  Pier  Horizontal  Stiffness') 
write (46,4780)  kvk.khk 

format (lx, 3x, fll . 1 , lx, 'kips/in ’, llx, fll . 1 , lx, 'klps/ln ’ ) 
write  (46,4800) 

format (lx,//, 20x, *  *  System  Parameters  and  Inputs  *’) 
write  (46,4850) 

format (lx, /, lx, ' Input  Forcing  Function  is  Horizontal  Component', 
+  '  of  the  1946  El  Centro’) 
write (46, 4875) 

format (lx, 20x, '  Earthquake  Acceleration  Time  History.’) 


write (46,4995) 

format (lx,/, lx, 'Vertical/Horizontal  Ground  Acceleration  Ratio' 
,3x,'Data  Time  Increment') 
write (46 , 4990)  amp,dtau 
format (lx. lOx, f6 . 3, tS5, f6 . 3, IX, ' sec ' ) 
write  (46,4900) 

format (lx,/, lx, 'Gravitational  Constant ', 3x, 'Percent  System 
Damping* ) 

write  (46,4950)  gravity , zeta* 100 . 

format (lx, 7x, f6 . 2 , lx, ' in/sec2 ' , lOx, f6 . 2 , lx, '%' ) 

write (46, 5000) 

format (lx,/, 25x, 'Mass  Matrix',/) 
do  5100  1=1,3 

write  (46,5050)  m(l.l)  ,m(i,2)  ,ro(i,3) 
format (lx, fl5.4, 5x. fl5. 4, Sx, fl5.4) 
continue 
write  (46,5200' 

format (lx,/, 25x, ’Damping  Matrix’ ,/) 
do  5300  1=1,3 

write  (46,5250)  cx(i.l)  ,cx(i,2)  ,cx(l,3) 
format  (Ix,fl5.4,5x,fl5.4,5x,fl5.4) 
continue 
write  (46 , 5400) 

format (lx,/, 25x, ’Stiffness  Matrix’,/) 
do  5500  1=1,3 

write  (46, 5450)  k(i,l),k(i,2),k(l,3) 

format (Ix,fl5.4,5x,fl5.4,5x,fl5.4) 

continue 

write (46 , 5700) 

format (lx,//) 

WRITE (46,6000) 

FORMAT (IX, 'Undamped  Natural  Frequencies' ,t3S, 'Mode  #l',tS0, 

'Mode  #2' ,t65, 'Mode  #3') 
write  (46,6001)  wl,w2.w3 

format (lx, t31 , f7 . 3, lx, 'rad/sec ' , t46 , f7 . 3, lx, 'rad/sec ’ , t62 , f7. 3, 

'  rad/sec') 

WRITE (46 . 6002) 

FORMAT (IX, 'Damped  Natural  Frequencies ', t35, 'Mode  #l',t50, 

'Mode  #2' ,t65, 'Mode  #3') 

WRITE (46,6500)  wl*sqrt (l-zeta**2) ,w2*sqrt (l-zeta**2) ,w3*sqrt (1-zet? 
format (Ix,t31,f7.3,lx. 'rad/sec ' , t46, f7 . 3, lx, 'rad/sec ' , t62, f7 . 3, 

'  rad/sec') 


endl  f 


319 


320 


11000 


20000 


If  (flag7.eq.l)  then 
f lagl0=f laglO+l 

write  (46,319)  x7 , y7 , t7, tlme7 

format  (lx, 'Side  block  crushing'  , t25, f9 .6, t35, f9 .6 , t50 , f9 .6 
+t6S,fS.2) 
endif 

If  (flag8.eq.l)  then 
flagl0=flagl0+l 

write  (46,320)  x8,y8,t8,time8 

format  (lx, 'Keel  block  crushing'  , t25, f9 .6, t3S, f9 .6, t50 , f9 . 6 
+  t6S, fS . 2) 
endif 

if (flaglO.le.O)  then 
write (46, 11000) 

format (lx, /, lx , 'No  failures  occurred.’) 
if (counter .eq. 1 .0  .and.  flaglO.le.O)  then 
go  to  60000 
endl  f 

1 f (counter -eq . 0 .0)  then 
ampacmax=ampacc 
ampacc=ampacc+ , 1 
counter=l . 0 

write (6,*)  'I  am  in  the  secondary  looping  stage.' 
endl  f 
endl  f 

i f  (ampacc . le . ampacmax)  go  to  20000 

if (counter .eq. 1 .0)  then 

ampacc=ampacc- . 0 1 

else  if  (counter .eq.O .0)  then 

ampacc=ampacc- . 1 

endif 

go  to  10000 
continue 


Since  the  vertical  equation  of  motion  is  uncoupled  from 
the  other  equations,  the  three  degree  of  freedom  system  need 
only  to  be  analyzed  to  obtain  damping  coefficients  for  both 
degree  of  freedom  mathematical  models.  Also,  the  maximum 
response  of  the  two  systems  will  be  determined  using  response 
spectrum  analysis  with  participation  factors  described  in 
Section  6.4. 

The  three  degree  of  freedom  equations  of  motions, 
undamped,  as  shown  in  Equations  4.4.11a,  b  and  c  are 
m^i  x  +  m13  c  +  k^  x  =  Xg 

m22  y  +  k22  y  =  ~m22  Yg 

and 

m33  -  +  m31  x  +  k3 3  -  =  -m31  xg 

where  m^  =  m22  =  M 

m13  =  m31  =  MKG 
m33  =  1^ 

kll  =  2Ksh  +  Kkh 
k22  =  2Ksv  +  Kkv 
k33  =  (B2/2 )  Ksv  -  W  KG. 

To  perform  modal  analysis,  consider  the  free  vibration  system 
m^  x  +  m13  6  +  kn  x  =  0 

m22  y  +  k22  y  =  0 

and 

m33  +  m31  x  +  k33  P  =  0. 

Assume  the  system  response  is  in  the  form 


x  =  (}>1  sin(wt  +  P  ) 
y  =  <j>2  sin  (cot  +  p ) 
and  8  =  <p3  sin(wt  +p )  . 


Now,  the  equations  of  motion  are 


mll  (“4>  l^2  sin(ajt  +  p  ))  +  m13  ( -tj>  ^2sin(o)  t  +  p)) 


+  k]^  (<J>  i  sin  (u>  t  +  p  ))  =  0 


m22  (“^ 2U ^  sin(wt  +P))  +  k22  (4>  2  sin (wt  +p))  =  0 


ra33  (“<f>  3  OJ2  sin(w  t  +  p  ))  +  m31  (-^  cu2  sin  (tut  +p)) 

+  k33  (<|>  3  sin  (to  t  +  p  ))  =  0. 


The  trivial  solution  to  the  above  set  of  equations  is 


sin  (out  +  p)  =0.  Assuming  sin  (cut  +p)  =  0,  then 


2  2 
kll-mllw  0  ~mi3  CO 


0  k22~m22W  0 


-m31(u 


0  k33  m33w 


(A. 3.1) 


These  equations  are  solved  for  w  by  setting  the  determinant  of 


the  first  matrix  equal  to  0.  However,  the  second  row  of  the 


matrix  is  uncoupled.  Thus  one  solution  for  is 


=  k22/m22, 


Now,  the  determinant  equation  reduces  to 


2  /mllkll+m33k33> 


kllk33 


CU  CO 


lmllm33  m3lm31, 


mllm33-m13m31 


Using  the  quadratic  formula,  the  other  values  of  u 2  are 


2  -bVb  -4c 
co  =  - 

2 


(A. 3. 2) 
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where  b  = 


-  k  "  y  ■ 


-  -  -  .  %  -  .  .  *.  \  \  \  %  v 


/'mllkll+m33k33\ 

\mllm33“ml3in3l/ 


and 


kllk33 

mllm33-ml3m31 


Using  parameters  found  in  Section  5  for  the  eleven  representa¬ 
tive  vessel-drydock  systems,  the  natural  frequencies  w  i,  o)2, 
and  u>3  for  each  system  is  calculated  and  listed  in  Table  A3.1. 
The  mode  shapes  are  determined  with  the  relationship 

(kn  -  mxl  u2)  -  m13u.'2  4>  3  =  0 

2 


or 


m 


13 


^ 3  kll  mH^2 


u  =  ,  ui  ^  • 


With  mode  shapes  determined,  the  equation  or  motion  can  be 
rewritten,  including  coupled  damping 

mu  *  +  mi3  §  +  cn  x  -t-  c13  e  +  kxl  x  -  0 

m22  y  +  c22  ir  +  k22  y  ■  0 

and 

..  • 

m3 3  0+  m33_  x  +  c33  6  +  c31  x  +  k33  =  0. 

At  natural  frequencies  w  2  andw3  the  coupled  three  degree 
of  freedom  system  acts  as  a  single  degree  of  freedom  system 
with  the  following  relationship: 

(4)  1/4)3)  x  =  9 

(01/413)  x  =  6 
and  (^1/^3)  x  =  6 . 

Using  these  relationships,  the  three  degree  of  freedom 
equations  of  motion  at  the  three  system  natural  frequencies 


TT 


*T> 


[mu  +  (^3/4>i)mi3]  x  +  [Cn  +  (<f>  3/<f>  i)  C31]  x  +  kn  x  =  0 

m22  y  +  c22  y  +  k22  y  ■  0 

and 

[  (^1/4)3 )m3i  +  m33]6  +  [  (4>  i/<t>  3)  c3i  +  c33]0  +  k33e  =  0. 
The  damping  coefficients  for  the  problem  are  5%  of  the 
critical  damping  for  each  mode, 

ten  +  (4*  3/4>  i)  C13 ]  =  2  CMi  wn  ,  wn  =  1,3  (A3. 3a) 

c22  =  2£  m22  “2  (A3. 3b) 

and 


[(4>l/4>3)^  c3i  +  c33]  =  2  M3wn,  un  -  1,3  (A3. 3c) 

where  ;;  =  percentage  of  critical  damping,  0.05 

M1  =  mll  +  ('ts/^iJn  m13 
and  M3  =  m33  +  (<t>  i/cj)  3 )  n  m3i. 

Substituting  known  values  and  rearranging  Equations  A3. 3a  and 
A3. 3c  yields  the  following  four  equations: 

Mode  1,  Cn  +  (<f3/4>i)i  ci3  =  2  Uia>i 

(4  i/4>  3)1  c3i  +  c33  =  2  £  M3  w  1 
Mode  2,  Cn  +  (<t>3/<t>i)3  ci3  =  2  C  Mi  w  3 

(<t>  3./$  3 )  3  c3 1  +  c3  3  ■  2  £  M3  w  3 
or  cn  +  (4  3/4>  1)  1  Ci3  =  2  C  Mi  w  1 

C11  +  (4*  3/^  l)  3  c13  *  2  C  Mj  u  3 
and  (4  1/4  3)i  c3i  +  c33  «  2  ^  M3  w  1 

(<t>  l/4>  3)3  C31  +  c33  “  2  £  M3  w  3 

Using  the  above  set  of  equations  and  equation  A3. 3b,  damping 
coefficients  Cn,  ci3,  c22,  c31'  and  c33  are  determined  and 
listed  in  Table  A3. 2  for  the  eleven  vessel-drydock  configura¬ 
tions.  Note  that  ci3  equals  to  c3i.  In  the  two  degree  of 
freedom  case,  the  damping  coefficients  cn,  ci2,  c2i,  and  c22 


>v>'. 

.  V  .>*, 


and  the  maximum  response  for  given  mode  is 


(6  ) 
max  n 


or  (6  ) 

max  n 


else  (6  ) 

max  n 


m 


11.62 


31 


[m33+(  3 )  nm3J 


for  to  co 


26.03 

(0 


m 


31 


[m33+(^1/<i,3)nm31] 


for  2.24 


331.53 


m31 

- -  for  to  , 

[m33+  (4>1/4>3)nm31] 


3  <  2.24 

<  a'l'  U3  <  12‘  74 
<*>3  >  12.74 


Since  being  uncoupled,  the  vertical  equation  of  motion  is 

»22  y  +  c2  t  +  k22  y  ■  -m22  yg 

c2y+k22y 


or 


y  + 


m 


22 


-y  =  -AMP  X 

g  g 


where  AMP  =  ratio  of  vertical/horizontal  El  Centro  earthquake 
accelerations . 

The  maximum  vertical  response  is 

(ymax)2  =  H.62  *  AMP  for  w2  <  2/24 

or 

(Ymax)2  =  26-‘ 03  *  AMP  for  2.24  <  « 2  <  12.74 

to 

else 

(ymax)2  -131^53  *  for  ^  >  12.74. 

to 

Now,  the  maximum  response  for  each  mode  of  the  x  and  6 
equations  and  for  the  uncoupled  y  equation  are  calculated. 

The  configuration  of  the  three  degree  of  freedom  system 
at  any  time  is  a  superposition  of  the  two  coupled  natural  mode 
shapes  along  with  the  independent  vertical  maximum  response. 
The  absolute  maximum  response  to  a  given  earthquake  is  the 
numerical  summation  of  the  maximum  response  of  each  mode  shape 
times  its  respective  participation  factor.  The  general 


formula  for  each  modal  participation  factor  is 
F  n  =  £  Mr  rn/  Mr  4>  2rn 

where  r  n  =  modal  particpation  factor  for  the  nth  mode 
r-  refers  to  which  equation  of  motion  is  being 
considered.  In  the  three  degree  of  freeddom 
case,  r  =  1  or  3 . 

Using  these  particpation  factors,  the  maximum  response  of  the 
three  degree  of  freedom  system  should  be  no  greater  than 

xmax  =  (xmax)l  ri  +  (xmax)3r  3' 
ymax  =  (ymax)2 

and  6  max  ■  ( emax)  1  rl  +  (6max)3r  3- 

These  maximum  response  are  determined  and  listed  in  Table  A3. 3 
for  the  eleven  vessel -drydock  configurations.  Once  again  due 
to  the  uncoupled  vertical  response,  the  maximum  response  of  x 
and  e  in  the  two  degree  of  freedom  is  identical  to  that  of  the 
three  degree  of  freedom  case. 


TABLE  A3.1 


Three  Degree  of  Freedom  Vessel-Drydock  System 
Natural  Frequencies 
(rad/ sec) 


SYSTEM 

W1 

co2 

w  3 

1 

5.454 

39.071 

37.230 

2 

3.883 

36.823 

34.367 

3 

4.013 

29.468 

23.791 

4 

2.819 

27.773 

21.511 

5 

2.848 

35.737 

32.742 

6 

9.690 

45.450 

40.863 

7 

7.953 

41.322 

36.351 

8 

6.976 

39.247 

34.045 

9 

5.863 

32.829 

29.060 

10 

6.857 

34.832 

29.921 

11 

6.177 

33.293 

28.010 

NOTE:  Two  Degree  of  Freedom  System  Natural  Frequencies 

correspond  to  and  oj3. 


TABLE  A3. 2 


Three  Degree  of  Freedom  Vessel-Drydock  System 
Damping  Coefficients 
(kim-sec/inch) 


SYSTEM 

C11 

c22 

c33 

c13 '  c3 1 

1 

65.398 

166.826 

1322882.125 

3920.416 

2 

58.612 

157.226 

938583.250 

2874.701 

3 

42.452 

125.824 

976657.438 

2829.964 

4 

37.308 

118.583 

682485.813 

2053.786 

5 

54.303 

152.588 

687504.500 

2159.187 

6 

188.657 

445.695 

7879946.000 

17130.227 

7 

167.100 

405.212 

6430934.500 

14269.135 

8 

155.786 

384.866 

5622059.000 

12660.506 

9 

43.963 

116.475 

1196244.250 

3340.509 

10 

33.833 

86.437 

836179.438 

2408.758 

11 

31.639 

82.618 

751062.813 

2185.568 

NOTE:  Two  Degree  of  Freedom  System  Damping  Coefficients 

correspond  to  cllf  c13,  c31  and  c33. 


0.33275 

0.10859 

0.01646 

0.27141 

0.12225 

0.02348 

0.65142 

0.19089 

0.02350 

0.53788 

0.21491 

0.02350 

0.21468 

0.12980 

0.03202 

0.45271 

0.08025 

0.00683 

0.46833 

0.09708 

0.00866 

0.46597 

0.10762 

0.01009 

0.64366 

0.15381 

0.01525 

0.65864 

0.13662 

0.01342 

0.67467 

0.14955 

0.01522 

Two  Degree  of  Freedom  System  Maximum  Response 
correspond  to  x  and  6. 


